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ABSTRACT 

This  report  describes  a  method  of  deriving  the  detection  performance  of  a  multistatic  sonar 
field  from  the  performance  of  a  field  of  similar  sonars  operated  monostatically.  It  allows  a 
direct  comparison  of  the  two  modes  of  operation,  thereby  quantifying  the  advantage,  if  any, 
of  multistatics.  The  method  is  derived  from  the  sonar  equation  in  the  noise-limited  regime. 
We  also  compare  two  different  network  architectures  for  tracking:  each  receiver  performing 
its  own  tracking  or  detection  information  being  passed  to  a  central  tracking  node.  We  start 
with  three  schematic  monostatic  detection-probability  curves,  ranging  from  almost  a  cookie 
cutter  to  an  exponential  shape,  which  has  a  long  low-probability  tail.  We  find  that  net¬ 
working,  whether  to  perform  multistatics  or  to  centralise  tracking,  brings  no  advantage  for  a 
cookie-cutter  detection  probability.  With  the  exponential  shape,  on  the  other  hand,  a  multi¬ 
static  field  can  be  spaced  at  about  twice  the  separation  for  the  same  detection  performance  as 
a  field  of  similar  sonars  operated  monostatically.  That  is,  a  given  area  can  be  covered  with 
about  one  quarter  the  number  of  sensors.  Centralising  the  tracking  allows  about  an  additional 
third  increase  in  sonar  separation. 
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Simple  Detection-Performance  Analysis  of 
Multistatic  Sonar  for  Anti-Submarine  Warfare 

Executive  Summary 

Multistatic  active  sonar  has  long  been  promoted  as  a  route  to  improved  performance  in 
anti-submarine  warfare.  However,  it  has  proven  difficult  to  quantify  the  benefits  of 
multistatics  in  operational  terms.  When  is  it  useful  and  when  is  it  not  worth  the  effort? 
This  report  introduces  a  simple  high-level  analysis  that  tackles  this  question  directly. 

We  work  from  the  point  of  view  that  operators  are  familiar  with  monostatic  sonar;  they 
know  what  to  expect  by  way  of  detection  performance  in  that  mode  of  operation.  The 
method  proposed  here  translates  monostatic  detection  range  into  multistatic  detection 
performance.  From  this,  we  can  obtain  the  separation  between  multistatic  sensors  that 
gives  the  same  detection  performance  as  similar  monostatic  sensors. 

The  simplicity  of  the  method  is  achieved  only  by  making  many  assumptions,  such  as, 
for  example,  the  availability  of  as  much  communications  capacity  as  may  be  required. 
Our  attitude  to  any  such  assumption  is  to  suppose  that  all  technical  obstacles  have 
been  overcome  and  to  ask:  what  then  is  the  tactical  benefit  of  multistatics,  in  quanti¬ 
tative  terms?  If  it  is  potentially  substantial,  then  this  may  provide  a  motivation  to  work 
on  the  technical  difficulties. 

To  give  a  direct  comparison  between  monostatics  and  multistatics,  we  adopt  a  field 
layout  of  collocated  source-receiver  pairs  arranged  on  a  triangular  grid.  We  consider 
large  sensor  fields  only;  conclusions  may  be  different  if  a  search  area  is  small  enough  to 
be  covered  with  a  few  sensors. 

Detection  range  is  measured  by  track-initiation  probability  rather  than  detection  pro¬ 
bability.  This  allows  us  to  compare  monostatics  and  multistatics  on  an  equal  footing,  as 
well  as  allowing  comparison  of  different  networking  architectures.  We  consider  the 
two  simplest  architectures:  each  receiver  performing  its  own  tracking  and  receivers 
passing  detection  information  to  a  central  tracking  node. 

The  results  show  that  networking,  whether  to  centralise  tracking  or  to  perform  multi¬ 
statics,  gives  very  little  advantage  when  the  monostatic  p ^  curve  is  close  to  a  cookie 
cutter.  This  has  long  been  understood,  but  a  quantitative  demonstration  has  not  been 
presented  before,  to  our  knowledge.  On  the  other  hand,  a  p curve  with  a  long  low- 
probability  tail  gives  the  opportunity  for  significant  multistatic  advantage:  sonobuoys 
can  be  placed  at  roughly  twice  the  separation  if  operated  multistatically  compared  with 
monostatic  operation.  That  is,  a  given  search  area  could  be  covered  with  about  one 
quarter  the  number  of  buoys.  This  applies  to  both  types  of  tracking  architecture,  with 
centralised  tracking  giving  about  a  further  third  increase  in  sonar  separation  over 
distributed  tracking. 

In  situations  where  the  false-track  rate  is  insupportably  high  for  fully  networked 
operation,  our  results  show  that  switching  to  distributed  tracking  brings  less  of  a 
performance  penalty  than  switching  to  monostatic  operation.  That  is,  when  multistatic 
centralised  cannot  be  done,  multistatic  distributed  is  better  than  monostatic  centralised. 
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1.  Introduction 


Anti-submarine  warfare  (ASW)  is,  and  always  has  been,  hard  (e.g.  [1,2]).  The  ADF  has  not 
given  it  priority  for  some  time  (e.g.  [3-5]): 

'Australian  ASW  capabilities  have  been  in  a  state  of  benign  neglect  for  a  couple  of 
decades  . . .  Other  considerations  took  precedence  over  ASW  when  setting  force  struc¬ 
ture  priorities.'  [3] 

Current  developments  in  submarine  technology  are  likely  to  make  ASW  harder  [6-8], 
while  the  on-going  proliferation  of  capable  submarines  in  our  region  (e.g.  [3,9])  suggests 
that  we  may  want  to  be  better  at  it.  This  situation,  coupled  with  the  nature  of  military 
leadership,  makes  it  natural  that  visions  for  the  way  forward  will  be  promulgated;  but 
how  to  assess  them?  There  is  a  dearth  of  actual  combat  experience  —  the  last  hot  ASW 
campaign  of  any  significant  length  took  place  during  World  War  2  [1,2]  — and  ASW  tech¬ 
nologies  are  complicated. 

Multistatic  active  sonar  is  an  example.  It  is  promoted  as  a  way  forward  for  ASW  and  it  is 
not  simple.  There  is  little  operational  experience  of  it  on  which  to  found  intuition  or 
military  judgement,  so  analysis  is  needed  to  help  assess  whether  it  is  indeed  a  way  for¬ 
ward  and,  if  so,  under  what  circumstances.  When  might  multistatics  be  useful  and  when  is 
it  not  worth  the  effort?  Answering  questions  like  these  is  hampered  by  an  'enduring 
problem  in  the  analysis  of  multistatic  systems  . . .  [namely]  the  lack  of  (even  approximate) 
closed-form  analytical  models  of  detection  performance'  [10].  This  report  introduces  a 
high-level  yet  quantitative  model  that,  though  not  closed-form,  aims  to  be  simple  enough 
to  be  readily  comprehensible  by  operators. 

The  concept  behind  the  analysis  involves  comparing  monostatic  and  multistatic  systems. 
Operators  are  familiar  with  the  behaviour  of  monostatic  active  sonar;  we  seek  to  translate 
this  experience  and  intuition  as  simply  and  directly  as  possible  to  a  'comparable' 
multistatic  system.  As  described  in  the  next  subsection,  the  study  question  is  formulated 
with  a  sonobuoy  field  in  mind,  but  the  analysis  applies  to  any  type  of  sonar  system.  We 
express  sonar-system  performance  in  terms  of  a  'range  of  the  day',  which  sweeps  up  all 
the  system,  environmental  and  target  details.  The  key  is  relating  the  range  achieved  by  a 
network  of  sonars  operated  multistatically  to  the  range  when  the  same  sonars  are  operated 
monostatically,  a  relationship  introduced  in  §1.2.2  and  detailed  in  §2.2.  The  success  of  the 
analysis  hangs  on  the  degree  to  which  this  relationship,  together  with  the  other  aspects  of 
the  method,  is  simple  enough  and  accurate  enough  for  the  purpose. 

The  level  of  simplicity  is  (and  can  only  be)  achieved  by  making  many  assumptions,  such 
as,  for  example,  the  availability  of  as  much  communications  capacity  as  may  be  required. 
It  is  also  assumed  that  methods  have  been  developed  to  handle  the  increased  false- 
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detection  rate  that  comes  with  networking  [11-14],  The  attitude  adopted  with  any  such 
required  assumption  is  to  suppose  that  all  technical  obstacles  have  been  overcome  and  to 
ask:  what  then  is  the  tactical  benefit  of  multistatics,  in  quantitative  terms?  If  it  is  found  to 
be  potentially  substantial,  then  this  may  provide  motivation  to  work  on  overcoming  the 
technical  difficulties. 

The  rest  of  this  section  elaborates  the  study  question,  sets  limits  on  the  scope  of  the  study 
and  explores  how  far  the  very  simplest  analysis,  using  cookie-cutter  detection  pro¬ 
babilities,  can  take  us.  This  sets  the  scene  for  the  analysis  proper,  the  method  of  which  is 
described  in  §2.  Section  3  gives  the  results  and  §4  contains  a  summary  and  a  statement  of 
conclusions.  Mathematical  details  are  largely  relegated  to  appendices,  which  also  contain 
the  MATLAB  code  used  to  generate  the  results. 

1.1  Study  Question  and  Scope 

The  following  are  cited  as  reasons  for  adopting  multistatics  (e.g.  [10,15-24]):  compared 
with  a  field  of  sonars  operated  monostatically,  the  same  field  operated  multistatically  may 

•  give  longer  detection  ranges 

•  increase  the  number  of  detection  opportunities  per  ping 

•  allow  higher  ping-repetition  rates 

•  complicate  the  tactical  situation  for  the  submarine 

•  where  more  than  one  receiver  makes  a  detection  simultaneously  (or  nearly  so), 
allow  improved  performance  in  localisation  (see  [25]  for  details),  classification  and 
tracking 

•  provide  a  tactically  significant  application  for  passive  sonobuoys  that  would  other¬ 
wise  be  of  only  marginally  utility. 

To  set  some  bounds  on  the  scope  of  the  work,  this  report  deals  only  with  the  first  two 
items,  which  are  linked  if  one  adopts  the  paradigm  that  defines  'detection  range'  in  terms 
of  the  range  dependence  of  a  cumulative  detection  (or  other)  probability  [26],  The  third 
and  last  points  are  treated  in  a  companion  report  [27];  the  fourth  and  fifth  points  are  left  to 
potential  future  work. 

In  formulating  the  study  question,  we  seek  to  draw  on  operator  experience  with  mono¬ 
static  sonar.  It  sets  a  standard;  operators  know  what  to  expect  by  way  of  detection 
performance  from  such  a  system.  Hence  we  frame  the  study  question  as: 

What  is  the  spacing  between  sonars  in  a  multistatic  field  that  gives  the  same  detec¬ 
tion  performance  as  a  monostatic  field  of  similar  sonars? 

To  make  this  work,  one  must  explain  how  to  determine  whether  two  sonar  fields  display 
the  'same  detection  performance'  and  what  is  meant  by  'similar  sonars'.  These  issues  are 
introduced  in  the  next  subsection  and  addressed  in  detail  in  §2.  Section  3  presents  results 
of  applying  the  method  to  a  field  of  collocated  source-receiver  pairs,  which  was  selected 
for  study  because,  of  all  possible  layouts,  it  allows  the  purest  comparison  between  mono¬ 
statics  and  multistatics.  Companion  reports  [27,28]  present  results  for  other  field  layouts. 
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1.2  Back-of-an-Envelope  Analysis  of  Sonar  Detection  Performance 

The  purpose  of  this  subsection  is  to  set  out  as  simply  as  possible  the  concept  on  which  the 
analysis  is  based.  We  begin  by  reiterating  the  simplest  model  of  detection  performance,  the 
'definite-range  law',  otherwise  known  as  the  'cookie-cutter'  detector.  This  provides  a 
pointer  to  a  simple  quantification  of  multistatic  detection  performance.  The  results  turn 
out  to  be  oversimplified,  for  reasons  explained  in  §1.2.3,  but  they  do  lead  to  some  general 
and  robust  conclusions  on  the  conditions  required  for  multistatics  to  be  beneficial. 


1.2.1  Monostatic  with  Definite-Range  Law 


The  'definite-range  law',  also  known  as  a  'cookie-cutter  detector',  is  the  simplest  model  of 
detection  performance.  In  its  most  basic  form,  the  'ideal  definite-range  law',(a)  detection 
probability  p$  as  a  function  of  range  R  to  the  target  is 


1 


Pd(R)  = 


if  R  <  Rq 
if  R  >  Rq  . 


(1) 


The  quantity  R0  *s  the  detection  range  or  'range  of  the  day'.  Equation  1  applies  to  an 
isotropic  detector,  that  is,  one  that  detects  equally  well  in  all  directions.  Many  types  of 
sonobuoy  are  good  examples. 


Figure  1  shows  three  identical  cookie-cutter  sensors  arranged  in  a  triangular  array.  As 
shown,  the  inter-sensor  spacing  Ds  is  not  optimum  because  there  is  a  region  covered  by  all 
three  sensors;  that  is,  the  overall  coverage  area  can  be  increased  without  introducing  a  gap 
by  choosing  a  somewhat  larger  Ds  value.  Simple  geometry  shows  that  the  maximum 
inter-sensor  spacing  for  which  there  is  no  gap  in  the  coverage  is 

Ds,max  =  ~  1-732R0.  (2) 

This  maximum  spacing  for  which  there  are  no  gaps  provides  a  simple  measure  of  the 
detection  performance  of  an  array  of  sensors. 


Figure  1:  Three  identical  sensors  positioned  at  the  vertices  of  an  equilateral  triangle  of  side  length 
Ds.  The  broken  lines  show  the  detection  ranges  R0  of  each.  The  maximum  spacing  for  no  gaps  in  the 
coverage  is  found  by  increasing  Dsfrom  the  value  shown  in  the  Figure  until  the  area  of  common 
overlap  of  the  three  detection-range  circles  just  goes  to  zero. 


(a) 


For  reference,  the  general  form  of  the  definite-range  law  is  given  in  Appendix  A.  1.1. 
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1.2.2  Relating  Monostatic  and  Multistatic  Ranges 

The  transmission-loss  term  of  the  sonar  equation  (e.g.  [29(ch.5)])  provides  a  way  of  con¬ 
necting  the  detection  ranges  of  a  field  of  sonars  operated  monostatically  with  the  same 
sonars  operated  multistatically.  In  the  monostatic  case  with  spherical  spreading  and 
neglecting  absorption,  the  transmission  loss  TL  is 

TLmono  =  40  k)g  R  (3) 

in  the  conventional  decibel  units.  The  factor  of  40  arises  because  the  acoustic  energy 
traverses  the  range  R  twice,  once  from  source  to  target  and  a  second  time  from  the  target 
back  to  the  receiver  that  is  collocated  with  the  source.  The  geometry  of  multistatic  detec¬ 
tions  is  shown  in  Figure  2.  Flere,  the  transmission  loss  is 

TLmulti  =  20  log  (Rst  Rtr)  .  (4) 

It  is  clear  that  the  two  transmission  losses  can  be  made  formally  equal  by  introducing  an 
equivalent  monostatic  range  Requiv: 

Requiv  =  \/%T  %R  •  (5) 

The  analysis  in  this  report  is  based  on  the  concept  of  using  Ret]uiv  to  reference  a  monostatic 
detection-probability  curve.  That  is,  if  a  monostatic  pd  curve  is  known  for  a  particular 
source,  collocated  receiver,  target  and  environment,  then  the  effect  of  separating  the  source 
and  receiver  while  keeping  everything  else  the  same  is  determined  by  using  Requiv  the 
monostatic  pd  curve.  This  analytical  concept  is  not  new;  it  is  implicit  in  many  studies  (e.g. 
[15,30])  and  has  also  been  explicitly  used  elsewhere  (e.g.  [10,24,31],  also  [32]  for  bistatic 
radar  applications).  The  basic  behaviour  of  the  model  is  illustrated  in  Figure  3  using  the 
definite  range  law  of  Equation  (1)  as  the  equivalent  monostatic  pd  curve.  The  resulting 
shapes  are  known  as  Cassinian  ovals  (e.g.  [33(§53)],  Appendix  A.2.1). 


Figure  2:  Geometry  of  multistatic  detections.  This  report  uses  only  the  distances,  not  angles  - 
sources  and  receivers  are  assumed  to  be  omnidirectional  and  the  aspect  dependence  of  the  target 
strength  is  neglected. 


4 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-TR-2562 


Figure  3:  Elementary  results  showing  the  geometry  of  Cassinian  ovals  -  detection  coverage  areas 
(green)  with  a  cookie-cutter  equivalent  monostatic  prf:  (a)  monostatic ,  (b-d)  bistatic  with  source- 
receiver  separation  shown.  Locations  of  the  receiver  are  shown  with  circles  and  the  source  with 
stars.  The  scale  bar  shozvs  the  length  corresponding  to  the  monostatic  detection  range  R0. 


The  method  allows  study  of  the  effect  of  'switching  on'  multistatics,  in  the  manner  implied 
in  the  study  question  of  §1.1.  It  can  be  applied  to  any  available  monostatic  pd  curve, 
whether  schematic  or  the  product  of  acoustic-propagation  calculations.  Because  it  is  a 
comparative  analysis,  the  impact  of  assumptions  is  diluted:  to  the  extent  that  they  affect 
the  monostatic  and  multistatic  situations  equally,  their  effects  cancel. 


1.2.3  Limitations  for  Treating  Multistatic  Networks  —  the  Way  Forward 

The  above  analysis  has  three  features  that  may  make  it  too  simple  to  be  useful.  These  are: 

•  use  of  the  ideal  definite-range  law, 

•  basing  coverage  area  on  single-ping  detection  curves,  and 

•  a  definition  of  Requiv  that  goes  to  zero  when  the  target  is  located  at  a  receiver. 

These  issues  are  discussed  in  the  following  paragraphs  and  ways  to  overcome  them  are 
developed,  thereby  establishing  a  basis  for  the  full  method  described  §2. 

Use  of  the  Ideal  Definite-Range  Law 

Suppose  the  sonar  field  comprises  collocated  pairs  of  sources  and  receivers,  and  receiver 
data  can  be  processed  either  monostatically  or  multistatically.  That  is,  there  is  both  a 
source  and  a  receiver  at  each  location  in  Figure  1,  and  we  can  choose  whether  each 
receiver  records  echoes  from  its  own  source's  pings  only,  or  listens  for  echoes  of  pings 
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originating  from  from  any  source.  If  the  monostatic  pd  curve  is  the  ideal  definite-range  law 
(Eqn  1),  then  it  is  clear  from  Figure  3  that  there  can  be  no  benefit  in  operating  the  field 
multistatically.  This  is  because  pd  is  unity  inside  the  detection  range  and  zero  outside,  so 
any  given  sensor  needs  no  assistance  for  targets  located  inside  its  detection  range,  and  can 
provide  no  assistance  outside  its  detection  range.  Benefit  from  multistatics  can  arise  only  if 
the  perfection  of  Equation  (1)  is  disrupted  in  some  way.(b)  Possibilities  are: 

•  Each  receiver  has  a  blind  zone  around  it,  arising  because  either  the  same  transducer 
is  used  both  to  transmit  and  receive,  which  cannot  be  done  simultaneously,  or  a 
receive  transducer  is  swamped  by  the  direct  blast  from  a  source  transducer. (c)  The 
blind  zone  is  a  region  of  zero  pL\  stretching  from  zero  range  out  to  a  value  related  to 
the  physical  length  of  the  ping  (details  in  §§2.1  and  2.2.2). 

•  No  realistic  pd  curve  is  unity  all  the  way  out  to  the  range  of  the  day.  Ranges  where 
pd  is  less  than  1.0  are  places  where  neighbouring  sensors  with  overlapping  detec¬ 
tion  regions  can  provide  mutual  support. 

•  The  abrupt  boundary  to  the  detection  region  is  also  unrealistic.  Real  pd  curves 
typically  tail  off  toward  zero  more  or  less  slowly,  and  usually  not  monotonically. 
The  tail  region  is  an  area  where  support  is  both  needed  from  and  can  be  supplied  to 
neighbouring  sensors. 

All  three  mechanisms  are  taken  into  account,  as  described  in  §2.1. 


Method  of  Defining  Coverage  Area 

Figures  1  and  3  suggest  an  analysis  with  a  measure  of  performance  related  to  detection 
coverage  area.  However,  both  figures  are  drawn  using  detection  ranges  derived  from 
single-ping  pd  values.  As  argued  elsewhere  [12,26],  it  is  better  to  base  a  definition  of 
detection  range  on  values  of  track-initiation  probability  pti.  (There  is  no  difference  between 
pd  and  pti  for  a  cookie-cutter  detector,  but  the  distinction  becomes  important  as  one  moves 
to  a  less  extremely  simplified  detection  model.) 

Track-initiation  probability  is  defined  in  terms  of  a  track-initiation  rule,  which  typically 
specifies  a  minimum  number  of  successful  detections  in  a  given  number  of  consecutive 
detection  opportunities.  The  concept  is  discussed  and  justified  in  §2.3,  and  details  sum¬ 
marised  in  Appendix  A.3.  We  use  the  3-in-5  track-initiation  rule  in  this  report.  The  full 
derivation  of  the  equations  is  given  in  a  companion  report  [34], 

In  summary,  defining  detection  range  through  pti  brings  two  main  benefits  over  defini¬ 
tions  in  terms  of  single-ping  pd  [26]: 

•  It  provides  a  level  of  recognition  of  the  impact  of  false  detections  arising  from 
oceanic  noise  [35].  ('False'  detections  of  non-target  objects  are  a  separate  issue.)  As 


(b)  These  arguments  apply  only  for  the  field  layout  considered  in  §3,  which  has  sources  and 
receivers  deployed  as  collocated  pairs.  The  effect  of  separating  sources  and  receivers  is  dealt  with  in 
companion  reports  [27,28];  it  can  give  multistatic  benefit  even  when  the  monostatic  pd  curve  is  an 
ideal  definite-range  law. 

(c)  The  term  'blast  zone'  is  also  used.  In  this  report,  we  use  'blind  zone'  to  cover  both  situations. 
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in  an  earlier  study  [12],  we  do  not  model  noise-induced  false  detections  explicitly, 
but  rather  rely  on  the  filtering  effect  of  the  track-initiation  step  to  mitigate  their 
effect. 

•  It  provides  a  means  of  quantifying  the  effect  of  networking  sensors.  This  is 
important  for  multistatics,  which  necessarily  involves  fields  of  geographically 
dispersed  sensors.  The  best  way  to  pool  detection  information  from  several  re¬ 
ceivers  is  an  important  question  that  can  be  addressed  in  terms  of  track  initiation; 
§2.3.2  explains  how. 


Definition  of  Requiv 

The  definition  of  Rec|uiv  in  Equation  (5),  though  consistent  with  the  sonar-equation  descrip¬ 
tion  of  spreading  loss,  has  the  undesirable  feature  of  going  to  zero  as  Rtr  — >  0.  It  means 
that  not  only  is  high  at  a  receiver  but  also  the  value  is  unaffected  by  the  source-receiver 
distance.  This  seems  inadequate;  an  improved  description  is  desirable  if  a  sufficiently 
simple  one  can  be  constructed.  The  following  are  possibilities: 

•  There  is  always  a  blind  zone  near  the  receiver,  which  means  that  p^  =  0  for  Rtr  near 
zero. 

•  Spreading  loss  changes  from  spherical  to  cylindrical  at  large  enough  range.  When 
the  source  is  far  from  the  receiver  but  the  target  is  close,  there  would  be  a  difference 
between  the  monostatic  and  multistatic  cases,  since  the  spreading  loss  would  be 
cylindrical  both  ways  for  monostatic  —  assuming  the  receiver  to  be  shifted  to  the 
source  location— but  cylindrical  one  way  (source  — »  target)  and  spherical  the  other 
(target  — »  receiver)  for  multistatic. 

•  Spherical  spreading  applies  in  the  far  field  of  radiated  sound.  Close  to  the  target, 
the  spreading  will  show  'near-field'  effects. 

•  Transmission  loss  is  caused  by  absorption  as  well  as  spreading.  The  absorption  loss 
involves  the  sum  of  RST  and  Rtr,  rather  than  its  product,  so  behaves  differently 
from  the  spreading  loss. 

We  examined  all  these  effects.  Near-field  and  cylindrical-spreading  effects  turn  out  to  be 
unimportant,  for  reasons  given  in  Appendix  B.l.  The  blind  zone  is  important  and  is 
included  in  the  results  of  §3,  but  it  is  an  unsatisfactory  resolution  of  the  Ret|uiv  problem 
because 

•  it  introduces  factors  unrelated  in  principle  to  the  heart  of  the  issue,  such  as  pulse 
length  and  whatever  level  of  pulse  compression  is  obtained  through  sonar-signal 
processing,  which  may  be  considerable  (§2.2.2),  and 

•  it  does  not  address  the  conceptual  problem  that  pL\  at  a  receiver  ought  to  show  some 
dependence  on  the  distance  to  the  source. 

Absorption  provides  the  answer,  as  illustrated  in  §2.2.1.  The  net  impact  on  the  results  of 
the  analysis  turns  out  to  be  small  for  the  geometry  and  distance  scales  of  interest  here,  so 
the  details  are  relegated  to  Appendices  B.1.1  and  B.1.2. 
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1.3  Summary  —  Interpreting  the  Study  Question 

The  analysis  of  §1.2  suggests  that  the  study  question  should  be  interpreted  in  the  following 
terms: 

•  'Similar  sonars'  means  that  the  performance  of  the  multistatic  field  is  determined 
by  using  Rec|uiv  as  the  range  in  a  monostatic  detection-probability  curve.  A  correc¬ 
tion  for  absorption  can  be  applied  if  desired. 

•  'Sonar  spacing'  means  the  maximum  spacing  for  which  the  area  covered  is  without 
holes  or  gaps. 

•  The  area  covered  should  be  defined  in  terms  of  track-initiation  probability  pti, 
rather  than  detection  probability  pd. 

The  last  two  points  indicate  that  'coverage  area'  as  defined  here  is  different  from  normal 
usage.  To  highlight  this,  we  give  it  a  new  name:  CATING,  for  'Coverage  Area  based  on 
Track  Initiation  and  No  Gaps'. 


2.  Method 


The  method  of  analysis  is  an  extension  of  earlier  work  [12]  to  incorporate  multistatics.  The 
measure  of  performance  is  the  maximum  sonar  separation  consistent  with  a  coverage  area 
without  gaps,  where  coverage  area  is  defined  in  terms  of  track-initiation  probability  pti 
(CATING).  The  use  of  pti  allows  us  to  quantify  the  effects  of  distributed  and  centralised 
tracking  in  a  manner  consistent  with  the  treatment  of  multistatics,  so  that  the  study  ques¬ 
tion  can  be  addressed  for  both  types  of  tracking.  This  is  explained  in  §2.3.2.  Before  this, 
§2.1  sets  out  two  simple  generalisations  of  the  definite  range-law  and  §2.2  details  how  this 
is  translated  to  multistatic  sensors.  Section  2.4  gives  the  method  of  calculating  the  measure 
of  performance. 


2.1  Models  of  Monostatic  Detection 


Although  the  method  of  relating  monostatic  and  multistatic  detection  performance  can  be 
used  with  any  type  of  monostatic  single-ping  pd  curve,  in  this  report  we  use  schematic 
curves,  in  part  to  keep  the  results  unclassified,  but  also  in  an  attempt  to  cover  a  reasonable 
range  of  detection  behaviour  with  a  few  curves. 


Two  classes  of  functions  are  used.  The  first  is  of  the  'Fermi  function'  type: 

1 


Pd(R)  = 


if  R  >  Rv 


1  +  1 0(K/K°-l)/& 

0  if  R  <  Rb , 


where 

•  R0  is  the  range  of  the  day,  defined  as  the  range  at  which  pd  =  0.5, 


(6) 
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•  b  is  a  'diffusivity'  parameter  that  describes  the  'tailing';  that  is,  rapidity  with  which 
probability  values  switch  from  near  unity  for  small  R  to  near  zero  for  large  R,  and 

•  Rb  is  the  radius  of  the  blind  zone. 

This  function  shows  quite  a  sharp  decrease  in  pd  at  ranges  around  R0,  at  a  rate  that 
depends  on  b.  As  b  — >  0,  it  approaches  the  definite-range  law  of  Equation  (1).  To  explore 
the  influence  of  the  diffusivity  parameter  on  sonar  field  performance,  we  performed  the 
analysis  twice,  using  the  two  values  of  b  shown  in  Figure  4. 


The  second  class  of  function  used  to  model  pd  behaviour  is  of  a  simple  exponential  type: 


Pd(R)  = 


'  10-0.30103R/R0 

0 


if  R  >  Rh 
if  R  <  Rb , 


(7) 


where  the  numerical  factor  is  required  to  ensure  that  pd(Ro)  =  0.5.  This  function  has  a 
much  longer  'tail  region'  where  0.1  <  pd  <  0.5  than  either  of  the  Fermi  functions  used,  as 
Figure  4  illustrates.  On  the  other  hand,  it  has  lower  pd  at  ranges  less  than  R0. 


The  length  Rb  of  the  blind  zone  is  related  to  the  physical  length  of  the  transmitted  acoustic 
pulse.  If  the  same  transducers  are  used  both  to  transmit  and  to  receive,  or  if  the  transmit¬ 
ted  waveform  does  not  allow  the  application  of  pulse-compression  techniques,  then  echos 
cannot  be  received  until  transmission  is  complete.  It  means  that  the  closest  detectable 
object  lies  at  a  range  equal  to  half  the  pulse  length.  Figure  4  is  drawn  with  Rb/ R0  =  0.075, 
which  corresponds  to  a  pulse  length  of  1.5  km  (pulse  duration  ~1  s)  if  the  range  of  the  day 
is  10  km. 


2.2  Multistatic  Detection  Probability 

It  is  a  foundation  of  this  analysis  that  the  multistatic  detection  probability  is  obtained  from 
either  Equation  (6)  or  (7),  as  the  case  may  be,  as  Pd(^equiv)'  where  Requiv  is  given  by  Equa¬ 
tion  (5).  There  are  two  provisos:  we  correct  the  p d  value  for  absorption  (though  the  effect  is 


Figure  4:  The  three  pd  curves  used  as  models  of  monostatic  detection  performance.  Range  is 
expressed  as  multiples  of  the  range  of  the  day  R0. 
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small),  and  the  treatment  of  the  blind  zone  is  different  in  the  multistatic  case;  it  is  not 
correct  simply  to  compare  Ret]mv  with  Rjr  These  two  issues  are  discussed  in  the  following 
subsections. 


2.2.1  Effect  of  the  Absorption  Correction 


As  described  in  Appendix  B.1.1,  the  absorption  correction  is  applied  to  a  signal  excess  (SE) 
rather  than  directly  to  a  pL\  value.  Hence,  the  method  of  correction  requires  converting  the 
value  of  Rd(^equiv)  to  an  SE,  applying  the  correction,  then  converting  back  to  give  the 
V d,corr  value. 


The  effect  of  the  procedure  is  displayed  in  Figure  5,  where  the  left  column  shows  un¬ 
corrected  values— i.e.  Rd(^equiv) —  and  the  right  column  shows  pd,corr-  In  this  figure  the 
monostatic  p^  values  are  from  the  Fermi  function  with  a  diffusivity  of  0.5,  the  attenuation 
coefficient  is  a  =  0.1  dB/km,  which  is  appropriate  for  an  acoustic  frequency  of  about  2  kHz 
[36],  and  the  range  of  the  day  is  set  to  10  km.  The  left  column  of  Figure  5  shows  the 
undesirable  behaviour  discussed  in  §1.2.3:  pd  rises  to  1.0(d)  at  the  location  of  the  receiver 
regardless  of  the  source-receiver  separation.  The  right-hand  panels  show  how  the  cor¬ 
rection  eliminates  the  effect:  pL\  values  peak  at  0.89  in  the  bottom  right  panel.  (This  is 
shown  more  clearly  in  Figure  6.)  It  is  a  rather  small  effect  in  view  of  the  large  range 
involved:  nine  times  the  range  of  the  day. 


uncorrected 


0  Rsr  =  9  R0  0 

corrected 


Figure  5:  Effect  of  the  correction  for  absorption:  left  column  -  uncorrected;  right  column  - 
corrected.  Contours  of  detection  probability  are  derived  from  a  Fermi- function  monostatic  p(^  of 
diffusivity  b  =  0.5  using  the  source-receiver  separations  shown,  a  range  of  the  day  of  10  km  and  an 
attenuation  coefficient  of  0.1  dB/km.  (The  contour  fill  colour  is  white  below  p(f  =  0.01.  Source  and 
receiver  locations  are  not  marked  in  the  manner  of  Figure  3  to  avoid  obscuring  prf  values  in  their 
immediate  vicinity.) 


(d)  The  peak  is  actually  0.9901,  because  that  is  the  value  of  the  Fermi  function  (Eqn  6)  at  R  =  0  when 
b  =  0.5.  (This  point  is  discussed  in  Appendix  A.1.2.) 
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Target  location  relative  to  receiver  (units  of  R0) 

Figure  6:  Detail  of  the  behaviour  of  near  the  source  or  receiver  location  in  the  two  bottom  panels 
of  Figure  5. 

2.2.2  Blind  Zone 

In  the  monostatic  case,  the  blind  zone  is  circular  with  a  radius  R b  equal  to  half  the  trans¬ 
mitted  pulse  length.  The  same  principle  applies  for  a  multistatic  geometry,  but  the 
separation  between  source  and  receiver  gives  the  blind  zone  an  elliptical  shape.  That  is, 
target  locations  such  that 

Rst  +  Rtr  <  Rsr  +  2Rb  (8) 

cannot  be  detected  because  their  echoes  arrive  at  the  receiver  at  the  same  time  as  some 
portion  of  the  pulse  directly  from  the  source.  Expression  (8)  with  the  inequality  replaced 
by  an  'equals'  sign  defines  an  ellipse  with  foci  at  the  source  and  receiver  and  the  shortest 
distance  from  either  equal  to  Rb.  Figure  7  shows  some  examples.  This  figure  contains  an 
important  message:  the  dramatic  increase  in  the  size  of  the  blind  zone  as  the  source- 
receiver  separation  increases.  Though  entirely  consistent  with  the  geometry  of  ellipses  (see 
Appendix  A.2.2),  it  nevertheless  appears  to  be  a  significant  negative  for  multistatics. 

The  large  size  of  the  piece  taken  out  of  the  covered  area  by  the  blind  zone,  so  apparent  in 
Figure  7,  provides  a  motivation  for  seeking  ways  around  it.  For  example,  it  is  possible  to 
design  a  pulse  waveform  so  that  careful  signal  processing  can  separate  an  echo  from  a 
partly  overlapping  direct  blast,  thereby  reducing  the  2 Rb  term  by  a  'pulse-compression 
factor'.  This  has  the  effect  of  increasing  the  eccentricity  of  the  ellipse,  which  means 
reducing  its  minor  axis  (Appendix  A.2.2).  In  other  words,  the  ellipse  shrinks  toward  the 
source-receiver  line.  The  effect  is  illustrated  in  Figure  8,  to  be  compared  with  Figure  7(c), 
which  shows  the  same  conditions  without  pulse  compression.  The  theoretical  maximum 
value  of  the  pulse  compression  factor  is  the  pulse  duration  multiplied  by  the  pulse  band¬ 
width  [15],  a  product  that  can  in  practice  exceed  several  hundred. 

The  benefits  of  pulse  compression  seem  clear,  but  are  in  fact  marginal  for  the  field 
geometry  studied  in  §3,  because  the  monostatic  capability  that  comes  with  collocated 
source-receiver  pairs  largely  covers  the  multistatic  blind  zones.  For  this  reason,  we  do  not 
consider  the  possibility  of  pulse  compression  in  the  rest  of  this  report. 
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Figure  7:  As  in  the  right  column  of  Figure  5  ,  but  with  the  inclusion  of  the  blind  zone,  which  is 
delineated  with  a  thin  black  line.  Locations  of  the  receiver  are  marked  with  circles  and  the  source 
with  stars,  except  in  panel  (a),  where  the  markers  would  obscure  the  blind  zone. 


(a)  pulse-compression  (b)  pulse-compression 


Figure  8:  Same  as  Figure  7(c),  but  with  two  values  of  the  pulse-compression  factor  (and  without  the 
node  markers). 


2.3  Coverage  Area  Based  on  Track-Initiation  Probability 

We  consider  a  track-initiation  rule  of  the  type:  start  a  track  if  there  are  at  least  m  detections 
in  n  consecutive  opportunities/6)  This  is  straightforward  for  a  single  monostatic  sonar,  as 


(e)  There  are  other  types  of  rule  (e.g.  [37,38]),  some  of  which  theoretically  give  a  lower  rate  of  false 
track  initiation.  However,  m-in-n  rules  remain  in  widespread  use  operationally. 
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outlined  in  §2.3.1,  but  questions  arise  with  networks  of  sensors.  If  there  are  many  sources, 
how  should  one  count  detection  opportunities  —  is  it  n  pings  from  one  source,  the  next  n 
pings  from  any  source,  or  n  pings  of  each  source?  We  adopt  the  last,  for  reasons  briefly 
summarised  in  §2.3.2  and  detailed  in  a  companion  report  [34],  When  the  field  has  many 
receivers,  does  each  perform  its  own  tracking,  or  is  detection  information  pooled  at  a 
central  tracking  node?  Should  the  answer  to  this  question  depend  on  whether  the  field  is 
operated  multistatically  or  not?  Exploration  of  these  two  questions  is  central  to  the  ana¬ 
lysis  in  this  report. 

2.3.1  Track-Initiation  Probability  for  a  Single  Sonar 

We  restrict  our  attention  to  the  3-in-5  track  initiation  rule.  That  is,  we  require  the  pro¬ 
bability  of  making  three  or  more  detections  in  five  consecutive  opportunities.  If  all  five 
opportunities  have  the  same  pd  value  and  are  statistically  independent,  then  the  track- 
initiation  probability  pd  is  [12] 

Fti  =  Fd3(10-15pd  +  6pd2).  (9) 

This  equation  is  plotted  in  Figure  9,  showing  that  pti  is  less  than  pd  when  pd  <  0.5  and 
greater  than  pd  for  large  pd  values.  The  effect  on  the  probability-range  curves  is  to  steepen 
the  transition  from  high  to  low  probability  values,  as  illustrated  in  Figure  10.  In  the  Fermi- 
function  case,  the  effect  of  using  pti  is  similar  to  using  a  pd  with  lower  diffusivity,  as  the 
blue  dotted  curve  in  Figure  10(b)  shows. 

The  pti  used  here  is  the  'local'  pti  of  reference  26.  It  is  preferred  to  a  cumulative  pti  because 
a  scenario  is  needed  to  calculate  cumulative  pti,  which  is  unnecessarily  restrictive  in  the 
present  context.  Also,  we  consider  the  3-in-5  rule  only;  some  insight  into  the  effect  of 
choosing  other  values  of  m  and  n  is  given  elsewhere  [12,34] . 


Figure  9:  Relationship  between  single-sensor  track-initiation  probability  and  detection  probability 
for  the  3-in-5  track-initiation  rule. 
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Range  (multiples  of  R0) 

Figure  10:  Comparison  of  the  range  dependence  of  detection  probability  (full  lines)  and  track- 
initiation  probability  (dashed  lines)  for  (a)  exponential  pd,  (b)  Fermi- function  with  b  =  0.5.  The 
blue  dotted  line  in  panel  (b)  shows  a  Fermi- function  with  b  =  0.25. 

2.3.2  Network  Architecture  and  Field  Track-Initiation  Probability 

In  applying  the  3-in-5  rule  to  a  networked  field  with  many  sources,  one  must  address  the 
question:  what  constitutes  5  detection  opportunities?  Our  view,  that  the  answer  should  be 
5  pings  of  every  source  in  the  field,  is  argued  in  detail  in  a  companion  report  [34].  In 
summary,  we  take  the  point  of  view  that  the  networking  ought  not  to  interfere  with  the 
detection  processing  of  the  sonars.  Each  sonar  processes  its  acoustic  data  in  the  same 
manner  whether  the  field  is  networked  or  not.  The  difference  lies  in  how  the  resulting 
detections  are  handled.  As  discussed  in  reference  34,  we  consider  four  types  or  'archi¬ 
tectures'  of  sonar  network: 

•  Monostatic  with  distributed  tracking:  The  sonars  operate  independently  of  each  other. 
Each  sonar  comprises  a  source-receiver  pair  and  each  performs  its  own  tracking 
based  solely  on  detections  from  its  own  pings.  A  given  sonar  must  ping  at  least  3 
and  up  to  5  times  before  it  knows  whether  to  start  a  track.  Any  data  fusion  that  may 
occur  consists  in  passing  information  on  initiated  tracks.  This  is  the  baseline  system. 

•  Monostatic  with  centralised  tracking:  Each  receiver  continues  to  respond  only  to 
echoes  from  its  own  source,  but  detection  information  is  passed  to  a  central  track¬ 
ing  node.  Each  source  still  pings  five  times,  so  that,  in  effect,  the  track-initiation  rule 
becomes  3-in-5J  for  a  field  comprising  /  sonars. 

•  Multistatic  with  centralised  tracking:  Each  receiver  in  the  field  processes  echoes  from 
any  source,  not  just  its  own,  with  detection  information  being  passed  to  a  central 
tracking  node.  Each  ping  produces,  in  principle,  /  detection  opportunities  for  a  field 
of  /  receivers.  If  there  are  also  /  sources,  each  of  which  pings  5  times,  then  the  track- 
initiation  rule  is  effectively  3-in-5/2.  Also,  if  there  are  3  or  more  receivers  in  the 
field,  then  it  is  possible  to  initiate  a  track  on  a  single  ping. 

•  Multistatic  with  distributed  tracking:  It  is  possible  that  tbe  false-alarm  rate  at  the 
central  tracking  node  in  the  centralised  multistatic  architecture  may  be 
insupportably  high.  In  an  attempt  to  mitigate  this,  it  has  been  suggested  that  a 
multistatic  field  could  return  to  distributed  tracking  [37,39],  That  is,  each  receiver 
operates  multistatically  but  does  its  own  tracking.  With  5  pings  from  each  of  / 
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sources,  this  returns  the  track-initiation  rule  to  3-in-5J.  Hence,  this  architecture 
should,  at  first  sight,  show  similar  sensitivity  to  the  false-detection  rate  as  a 
monostatic  field  with  centralised  tracking. 

To  provide  a  basis  for  quantifying  the  detection  performance  of  fields  of  sensors  in  a  way 
that  treats  different  architectures  on  an  equal  footing,  we  need  to  generalise  the  concept  of 
track-initiation  probability  in  §2.3.1  to  encompass  the  whole  field.  The  required  mathema¬ 
tical  derivations  are  detailed  in  reference  34,  with  results  listed  in  Appendix  A.3. 

Figure  11  displays  the  impact  of  network  architecture.  The  centralised  multistatic  case 
(panel  d)  has  field  track-initiation  probability  ptif  exceeding  0.75  over  the  whole  area 
shown  (more  than  0.85  over  most  of  the  area),  whereas  the  distributed  monostatic  case  in 
panel  (a)  shows  many  regions  with  pti  f  below  0.15. 

The  comparison  between  centralised  monostatic  and  distributed  multistatic  cases  in  panels 
(b)  and  (c)  is  interesting  in  view  of  the  argument  above  that  these  two  should  have  similar 
sensitivity  to  false  detections.  For  the  conditions  shown  in  Figure  11,  the  distributed  multi¬ 
static  architecture  clearly  performs  better,  in  the  sense  that  the  pti  f  values  between  the 
sensors  remain  above  0.4  everywhere  in  Figure  11(c),  but  fall  below  0.25  in  several  places 
—  the  red  patches  — in  Figure  11(b). 
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Figure  11:  Contours  of  field  track-initiation  probability  in  the  interior  of  a  large  sonar  field  for  the 
four  network  architectures.  The  field  comprises  posts  (collocated  source-receiver  pairs),  indicated  by 
the  small  dots  at  the  centre  of  each  monostatic  blind  zone,  on  a  triangular  grid  spaced  at  2.5  times 
the  range  of  the  day.  The  monostatic  p^  curve  is  a  Fermi  function  with  diffusivity  b  =  0.5  and  the 
absorption  correction  is  applied  to  the  multistatic  cases.  The  effect  of  the  multistatic  blind  zones  is 
discernable  in  panels  (c)  and  (d)  as  sharp-edged  structures  in  the  lower- ptij  regions  between  the 
sensors,  though  it  is  combined  with  the  range  dependence  of  Fti,f-  The  y>ti,f=  0.5  contour  is  de¬ 
lineated;  in  (d),  p fi.fis  nowhere  less  than  0.75  over  the  area  shown. 
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2.4  Measure  of  Performance 

In  view  of  the  study  question  (§1.1)  and  its  interpretation  (§1.3),  the  measure  of  perfor¬ 
mance  (MOP)  adopted  is  the  maximum  separation  between  sensors  such  that  pti j  is  not 
less  than  0.5  in  the  interior  of  a  large  field.  Thus,  the  sensors  in  Figure  11(d)  are  too  closely 
spaced;  those  in  all  the  other  panels  of  Figure  11  are  too  far  apart.  To  determine  the  MOP 
value  in  a  given  case,  the  inter-sensor  spacing  is  adjusted  manually  until  pti  f  =  0.5  contours 
just  appear  in  between  the  central  sensor  and  the  ring  of  sensors  around  it,  as  illustrated  in 
Figure  12.  We  carried  out  this  procedure  to  the  nearest  0.1  R0  for  the  exponential  pL\  curve 
and  the  nearest  0.01  R0  for  the  Fermi-f unction  cases  like  Figure  12. 

In  determining  an  MOP  value.  The  monostatic  blind  zone  is  excluded  from  the  require¬ 
ment  that  pti  f  be  not  less  than  0.5.  That  is,  a  small  region  of  zero  pti  f  surrounds  each 
sensor,  as  shown  in  Figure  12.  We  do  not  reduce  the  inter-sensor  spacing  so  far  as  would 
be  required  to  fill  this  in. 

As  emphasised  in  the  next  section,  the  analysis  is  intended  to  apply  to  large-area  search 
only.  The  implication  for  MOP  determination  is  that  we  must  use  fields  extensive  enough 
for  there  to  be  no  edge  effects  near  the  central  sensor.  For  example,  edge  effects  appear  in 
Figure  11(d)  as  the  very  small  pale-green  patches  in  the  corners  of  the  panel.  These  lie  well 
away  from  the  central  sensor.  No  edge  effects  are  discernable  in  Figures  ll(a-c).  In  each 
case  analysed,  we  checked  for  edge  effects  by  adding  successive  rings  of  sensors  until  ptij 
values  in  the  vicinity  of  the  central  sensor  were  unaltered  by  the  new  ring  of  sensors.  More 
rings  of  sensors  are  needed  for  the  exponential  monostatic  pL\  curve  than  the  cookie-cutter 
(low-diffusivity  Fermi),  as  could  be  expected. 


Figure  12:  Like  Figure  11(d),  but  showing  the  region  of  the  central  sensor  for  inter-sensor  sepa¬ 
rations  of  (a)  2.73  R0  and  (b)  2.74  R0.  The  small  black  dots  on  the  edges  of  the  yellow  areas  in  (b) 
show  regions  where  p#  * just  falls  below  0.5. 
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3.  Results 


For  convenience,  we  first  collect  all  the  assumptions  and  parameter  values  used,  then 
show  examples  of  the  patterns  of  pti  t  and  present  the  results. 

3.1  Summary  of  Modelling  Paradigm  and  Parameter  Values 

•  Study  question:  What  is  the  spacing  between  sonars  in  a  multistatic  field  that  gives 
the  same  detection  performance  as  a  monostatic  field  of  similar  sonars? 

•  Field  layout:  Triangular,  as  in  Figure  11,  with  a  post  (collocated  source-receiver  pair) 
at  each  sensor  location. 

•  Network  architectures  considered: 

♦  monostatic  distributed:  Each  receiver  responds  only  to  echoes  of  pings  from  its 
collocated  source,  and  each  receiver  performs  tracking  based  only  on  its  own 
detections.  This  is  the  baseline  case. 

♦  multistatic  distributed:  Receivers  register  echoes  of  pings  from  any  source,  but 
each  receiver  still  performs  its  own  tracking. 

♦  monostatic  centralised:  Each  receiver  responds  only  to  echoes  of  pings  from  its 
collocated  source.  Detection  information  is  passed  to  a  central  tracking  node. 

♦  multistatic  centralised:  Receivers  attempt  detections  on  every  ping  and  pass  the 
information  to  a  central  tracking  node. 

•  Models  of  monostatic  detection  probability  vs.  range  curves: 

♦  Fermi,  b  =  0.1,  ♦  Fermi,  b  =  0.5,  and  ♦  exponential. 

All  are  parametric  in  R0,  the  range  of  the  day,  so  this  need  not  be  specified. 

•  Multistatic  detection  probability:  From  monostatic  pd  curves  using  Requiv  =  V^ST%R  • 

•  Absorption  correction:  Applied  using: 

♦  attenuation  coefficient  of  0.01  dB  per  R0,  and 

♦  noise  standard  deviation  of  8  dB. 

•  Blind  zone:  Applied  using: 

♦  transmitted  pulse  length  of  1.5  km,  taken  to  be  0.15R0  (hence  =  0.075Ro),  and 

♦  no  pulse  compression. 

•  Measure  of  performance:  maximum  inter-sensor  separation  such  that  field  track-initi¬ 
ation  probability  in  the  interior  of  a  large  sensor  field  is  nowhere  less  than  0.5, 
except  for  the  residual  blind  zones  in  the  immediate  vicinity  of  each  sensor. 


3.2  Patterns  of  Field  Track-Initiation  Probability 

Figure  11  shows  patterns  of  pti  f  for  the  case  of  Fermi,  b  =  0.5.  Corresponding  results  for  the 
other  two  monostatic  pL\  curves  are  shown  in  Figures  13  and  14.  In  Figure  13,  the  mono¬ 
static  pd  curve  is  close  to  a  cookie  cutter,  so  the  networking  benefit  is  small,  as  suggested  in 
§1.2.3.  Figure  14  shows  the  opposite  extreme:  the  monostatic  pL\  curve  has  a  long  tail  where 
Pd  lies  in  the  range  ~0. 1-0.5,  which  is  the  best  situation  for  gaining  advantage  through 
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Figure  13:  Like  Figure  11,  but  for  the  Fermi  monostatic  curve  with  dijfusivity  b  =  0.1  and  an 
inter-sensor  separation  ofl.9R0 
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Figure  14:  Like  Figure  11,  but  for  the  exponential  monostatic  curve  and  an  inter-sensor  sepa¬ 
ration  of  6  R0 
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networking  [12],  Figure  14(d),  the  multistatic  centralised  case,  has  ptiy  over  0.95  every¬ 
where,  apart  from  in  the  residual  blind  zones,  so  the  sensors  can  be  spaced  much  further 
than  the  6  R0  shown.  Figure  15  illustrates  the  determination  of  the  maximum  spacing  for 
this  case. 


3.3  Multistatic  Advantage  for  a  Field  of  Collocated  Sources  and 
Receivers 

The  results  of  the  analysis  are  presented  in  Figure  16.  The  quantitative  performance  gain  in 
operating  the  field  multistatically  can  be  gauged  by  comparing  a  blue  bar  with  the  red  bar 
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Figure  15:  Like  Figure  14(d),  but  for  inter-sensor  separations  of  (a)  8.1  R0  and  (b)  8.2  R0.  The 
regions  with  pf^  below  0.5  in  panel  (a)  are  due  to  edge  effects. 

to  its  left.  The  inter-sonar  separation  of  V3  Rq  appropriate  for  a  cookie-cutter  pd  curve  is 
shown  for  comparison. 

The  left-most  group  of  bars  in  Figure  16  relate  to  a  monostatic  p ^  curve  that  is  almost  a 
cookie-cutter.  The  multistatic  advantage  is  very  small  with  distributed  tracking  (solid  bars) 
and  not  much  larger  for  centralised  tracking  (striped  bars). 

In  contrast,  the  right-most  group  of  bars,  for  the  exponential  monostatic  pL\  curve,  show 
significant  multistatic  advantage  for  both  tracking  architectures.  For  the  centralised 
architecture  (striped  bars),  the  maximum  inter-sensor  spacing  for  the  multistatic  field  is  a 
little  under  double  that  of  the  monostatic  field  (8.1  Rq  cf.  4.8  Rq)  and  for  distributed  track¬ 
ing  it  is  more  than  double,  at  6.3  Rq  compared  with  2.8  Rq.  As  to  the  impact  of  tracking 
architecture,  for  both  monostatic  and  multistatic  fields  centralised  tracking  does  about  a 
third  better  than  distributed  tracking.  In  summary,  using  the  monostatic  distributed 
configuration  as  a  baseline,  moving  to  multistatics  approximately  doubles  the  field  per¬ 
formance  metric  and  centralising  the  tracking  increases  it  by  about  another  third. 

The  middle  group  of  bars  show  a  case  of  a  monostatic  pL\  curve  that  is  clearly  not  a  cookie 
cutter,  but  with  a  tail  region  of  relatively  limited  extent.  Multistatic  advantage  is  present 
but  modest.  This  time,  again  starting  with  the  monostatic  distributed  baseline,  moving  to 


b  =  0.1  5  =  0.5 
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Figure  16:  Values  of  the  measure  of  performance  for  the  12  cases  calculated.  Each  pair  of  red  and 
blue  bars  shows  the  multistatic  advantage  for  that  combination  of  monostatic  p^  curve  and  tracking 
architecture.  The  baseline  network  configuration  is  monostatic  with  distributed  tracking,  shown  as 
the  left-most  bar  in  each  of  the  three  groups. 
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multistatics  increases  the  performance  metric  by  about  20%  and  centralising  the  tracking 
adds  something  over  another  10%. 

The  magnitude  of  the  tailing  in  the  monostatic  pL\  curve  can  be  quantified  by  'tail  width' 
PVt,  which  we  define  as  the  distance  from  the  range  at  which  pd  =  0.5  to  the  range  at  which 
pd  =  0.1  (Appendix  A.1.3).  Figure  17  shows  the  data  of  Figure  16  plotted  against  tail  width. 
The  extent  to  which  the  exponential  represents  extreme  behaviour  is  apparent.  It  is  also 
clear  that  multistatic  advantage,  which  is  the  difference  between  corresponding  red  and 
blue  curves,  grows  with  increasing  tail  width.  There  is  no  such  advantage,  nor  any  advan¬ 
tage  in  centralising  detections  for  the  purpose  of  tracking,  for  a  monostatic  pd  curve  with 
zero  tail  width,  such  as  a  cookie  cutter.  (f) 

A  separate  question  concerns  the  relative  performance  of  the  centralised  monostatic  and 
distributed  multistatic  architectures  (middle  two  bars  of  each  group  in  Figure  16  or  middle 
two  lines  in  Figure  17).  This  becomes  interesting  if,  for  some  reason,  one  cannot  cope  with 
the  false-track  rate  of  the  fully  networked  case.  In  such  a  situation,  is  it  better  to  give  up 
centralised  tracking  or  multistatics?  In  every  case  examined,  distributed  multistatics  out¬ 
performs  centralised  monostatics,  as  Figure  17  makes  clear,  though  the  difference  is 
marginal  if  the  tail  width  is  small.  Other  considerations,  such  as  the  high  sensitivity  of  a 
multistatic  field  to  sonar  positioning  errors,  may  cause  monostatics  to  be  favoured  in  some 
situations. 


Figure  17:  Same  as  Figure  16,  but  plotted  against  tail  width -the  distance  between  the  range  at 
which  pd  =  0.5  and  the  range  at  which  pd  =  0.1  (Appendix  A.1.3).  The  definite-range  law  ('cookie- 
cutter')  has  zero  tail  width.  All  curves  converge  to  a  maximum  inter-sensor  spacing  o/V3R0  (-1.73 
R0)  as  the  tail  width  goes  to  zero. 


(0  This  conclusion  applies  only  for  the  ideal  definite-range  law  of  Eqn  (1),  where  there  is  detection 
certainty  inside  the  detection  range.  If  pd  is  not  1.0  inside  the  cookie  cutter  (e.g.  as  in  Appendix 
A.1.1)  then  networking  will  bring  an  advantage  even  zero  tail  width.  Reference  12  contains  an 
example.  See  also  the  comments  in  footnote  (b):  multistatic  advantage  can  also  be  obtained  in  envi¬ 
ronments  with  zero  tail  width  by  separating  the  sources  and  receivers,  or  adding  extra  receivers  to  a 
field  of  posts. 
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4.  Summary,  Comments  and  Conclusion 

This  report  describes  a  modelling  method  for  comparing  the  detection  performance  of  a 
field  of  multistatic  sonar  systems  with  that  of  a  field  of  similar  sonars  operated  mono- 
statically.  The  basis  of  the  method  is  the  sonar-equation  description  of  transmission  loss 
due  to  spreading,  which  leads  to  a  very  simple  connection  between  the  multistatic  ranges 
and  an  equivalent  monostatic  range: 

Kequiv  =  VR STRTR  • 

It  means  that  multistatic  detection  probabilities  can  be  obtained  directly  from  monostatic 
pd  values.  The  modelling  paradigm  supposes  that  the  monostatic  pd  values  are  known  a 
priori,  so  the  comparison  is  straightforward.  For  simplicity  and  to  keep  the  study  un¬ 
classified,  we  use  three  simple  monostatic  pd  curves:  one  that  is  almost  a  cookie  cutter,  one 
with  a  long  tail  where  pd  lies  in  the  range  ~0. 1-0.5,  and  a  third  intermediate  case. 

The  form  of  the  sonar  equation  leading  to  the  above  definition  of  Requiv  applies  in  the 
noise-limited  regime.  Acoustic-propagation  calculations  for  the  reverberation-limited 
regime  show  as  much  regularity  as  when  noise-limited  (e.g.  Fig.  7  in  ref.  30),  so  an  analysis 
along  the  lines  of  the  analysis  presented  in  this  report  should  be  possible  by  starting  with 
the  reverberation-limited  form  of  the  sonar  equation  (e.g.  [11]).  Of  course,  the  details 
would,  and  conclusions  might,  be  different  from  those  presented  here. 

The  results  reported  in  §3.3  include  the  effect  of  the  blind  zone  between  source-receiver 
pairs  and  a  correction  for  the  effect  of  transmission  loss  due  to  absorption,  using  a  value  of 
the  attenuation  coefficient  appropriate  to  a  sonar  frequency  of  ~2  kFIz,  although  both  turn 
out  to  be  minor  given  the  chosen  field  layout.  Other  potential  corrections  —  near-field 
effects  and  the  effect  of  a  transition  from  spherical  to  cylindrical  spreading— were  investi¬ 
gated  and  found  to  be  negligible. 

To  give  the  most  direct  comparison  possible  between  monostatics  and  multistatics,  we 
adopted  a  field  layout  of  collocated  source-receiver  pairs  ('posts')  arranged  on  a  triangular 
grid  with  equal  spacing  between  each  post.  We  considered  large  sensor  fields  only;  con¬ 
clusions  may  be  different  for  a  search  area  small  enough  that  it  can  be  covered  with  a  few 
sensors.  For  example,  Travaglione  and  Forward  consider  a  single  source  surrounded  by  a 
circle  of  2-5  equally  spaced  receivers  [31].  They  find  an  optimum  source-receiver  spacing 
of  ~1.3Rq,  but  their  analysis  differs  from  ours  in  more  than  just  the  field  layout:  they  base 
their  coverage-area  metric  on  pd  rather  than  and  they  assume  very  short  acoustic  pulses 
(50  ms),  so  that  the  blind-zone  areas  are  relatively  small. 

Aspect-dependent  target  strength  could  be  included  in  the  present  analysis  in  a  similar 
manner  to  the  absorption  correction.  We  have  not  investigated  this.  Although  the  possi¬ 
bility  of  'glint'  detections  is  held  to  be  an  advantage  of  multistatics  over  monostatics,  the 
effect  is  probably  not  as  significant  when  sources  and  receivers  are  collocated  compared  to 
when  they  are  separated,  particularly  if  the  submarine  CO  is  aware  of  (or  has  a  shrewd 
guess  at)  the  layout;  for  then,  in  addition  to  knowing  the  locations  of  the  sources,  he  or  she 
also  knows  where  all  the  receivers  are. 
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Perhaps  the  most  significant  shortcoming  of  the  analytical  method  is  the  neglect  of  time,  a 
feature  that  is  typical  of  methods  based  on  coverage  area.  We  take  no  account  of  the  time 
required  for  five  pings  from  all  sources  in  the  field,  nor  of  submarine  motion  during  this 
time.  On  the  plus  side,  it  allows  us  to  avoid  the  difficult  issue  of  ping  sequencing  (which 
others  have  considered  [40,41]),  but  this  scarcely  compensates  for  the  negatives  that  the 
approximation  entails.  The  best  hope  for  mitigating  its  effects  may  be  future  technical 
developments  that  render  its  impact  less  significant,  such  as  methods  for  allowing  simul¬ 
taneous  pings  by  several  sources  (e.g.  [40]),  or  continuous  active  sonar  (e.g.  [42-44]),  which 
would  eliminate  the  need  for  ping  sequencing  altogether. 

The  measure  of  performance  is  the  maximum  inter-sensor  spacing  such  that  the  field 
track-initiation  probability  between  sensors  is  not  less  than  0.5.  The  use  of  track-initiation 
probability  fills  two  purposes: 

•  Since  track  initiation  acts  as  a  false-detection  filter,  use  of  pti  f  provides  some  level  of 
robustness  against  the  effect  of  false  detections,  which  otherwise  are  not  considered 
in  this  analytical  method. 

•  It  provides  a  way  of  quantitatively  comparing  the  performance  of  different  network 
architectures,  both  monostatic  versus  multistatic,  and  also  distributed  versus  cen¬ 
tralised  tracking. 

It  should  be  carefully  noted  that  the  measure  of  performance  relates  to  detection  performance 
only.  We  take  no  account  of  other  issues  in  multistatics,  such  as  localisation  error  or  the 
time  required  to  achieve  a  given  level  of  confidence  in  the  identity  of  the  contact. 

The  results  show  that  networking,  whether  to  centralise  tracking  or  to  perform  multi¬ 
statics,  gives  very  little  advantage  when  the  monostatic  pL\  curve  is  close  to  a  cookie  cutter. 
This  has  long  been  understood,  but  a  quantitative  demonstration  has  not  been  presented 
before,  to  our  knowledge.  On  the  other  hand,  the  long-tail  pc\  curve  gives  the  opportunity 
for  significant  multistatic  advantage:  with  the  exponential  curve,  sonobuoys  could  be 
placed  at  roughly  twice  the  separation  if  operated  multistatically  compared  with  mono¬ 
static  operation.  That  is,  a  given  search  area  could  be  covered  with  about  one  quarter  the 
number  of  buoys.  This  is  true  for  both  types  of  tracking  architecture,  with  centralised 
tracking  giving  about  a  third  better  performance  than  distributed  tracking  (Figs  16, 17). 

In  situations  where  the  false-track  rate  is  insupportably  high  for  fully  networked  opera¬ 
tion,  our  results  show  that  switching  to  distributed  tracking  brings  less  of  a  performance 
penalty  than  switching  to  monostatic  operation.  That  is,  when  multistatic  centralised  can¬ 
not  be  done,  multistatic  distributed  is  better  than  monostatic  centralised. 

In  terms  of  providing  general  conclusions,  the  field  layout  considered  here  is  quite  restric¬ 
tive-collocating  sources  and  receivers  is  perhaps  not  the  best  way  of  employing  a  multi¬ 
static  capability.  We  have  extended  the  analysis  to  other  field  layouts,  the  results  of  which 
appear  in  companion  reports  [27,28], 
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Appendix  A: 

Mathemaical  Details 

This  appendix  collects  some  results  in  support  of  the  analysis  in  the  body  of  the  report. 


A.l.  Variants  and  Properties  of  the  Detection-Probability  Models 

A.1.1  Definite-Range  Law  and  Exponential 


Equation  (1)  is  referred  to  as  the  'ideal  definite-range  law'.  The  general  version  is 


Pd(R)  = 


jpo 

jo 


R<R0 

R>R0. 


(At) 


The  only  effect  of  this  generalisation  is  that  R0  does  not  fulfil  our  definition  of  'range  of  the 
day'  if  p0  is  chosen  to  be  less  than  0.5.  In  such  a  case  there  is  no  range  of  the  day. 


In  reference  12,  the  exponential  model  is  written 

Pd(R)  =  PoexP(-R/«)-  (A2) 

For  p0  =  1,  this  gives  the  same  values  as  Equation  (7)  — the  blind  zone  aside  — if  a  is  set  to 


a  = 


Ro 

lnlO  log  2 


1.44271*0 


(A3) 


If  it  is  desired  to  generalise  Equation  (7)  by  introducing  a  p0  value,  then  one  should  write 

pd  (R)  =  p0 10“log(2p°  WR° ,  (A4) 


so  that  R0  remains  the  range  of  the  day  (i.e.  the  range  at  which  pd  =  0.5).  This  formulation 
requires  p0  >  0.5;  for  otherwise  there  is  no  range  of  the  day.  To  make  Equation  (A2)  give 
the  same  values  as  Equation  (A4),  set 


In 10  log(2p0)  ‘ 


(A5) 


A.1.2  Fermi  Function 


In  reference  12,  the  Fermi  model  is  written,  apart  from  some  minor  changes  of  notation  to 
suit  present  purposes,  as 


Pd(R)  =  R0 


l  +  exp(-r0/fl) 


1  +  exp 


R~r0 

V  a 


(A6) 
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(Note  that  P0  ¥=  p0  and  r0  +  R0,  as  explained  below.)  The  numerator  ensures  that  pd  =  P0  at 
R  =  0,  but  it  complicates  the  interpretation  of  the  r0  parameter;  for  the  range  of  the  day  R0 
is 


R0  =  r0  +  a  ln 


+  e_r°/fl 


(A7) 


This  cannot  be  inverted  to  obtain  an  expression  for  tq  in  terms  of  R0.  In  the  present  report, 
we  choose  to  drop  (l  +  e_'0//<7 )  factor  in  the  numerator.  It  means  that  pd( 0)  +  p0,  but  the 
value  of  pd(0)  not  accessible  anyway  because  the  region  around  R  =  0  is  covered  by  the 
blind  zone. 


Equation  (6)  can  be  generalised  to  include  a  p0  factor  while  retaining  the  interpretation  of 
R0  as  the  range  of  the  day  by  writing,  again  apart  from  the  blind  zone. 


Pd(R)  = 


_ Po _ 

l  +  l0{R[i+i’log(2Po-1)]-P°}/^  ' 


(A8) 


It  is  apparent  that  this  reduces  to  Equation  (6)  when  p0  =  1.  Equation  (A8)  applies  only  for 
p0  >  0.5;  for  otherwise  there  is  no  range  of  the  day.  To  make  Equations  (A6)  and  (A8)  give 
the  same  values,  set 


?0  l  +  &log(2p0-l)' 


a  = 


br0 
In  10  ' 


Po=- 


Po 


1  +  10 


-l/b  ‘ 


(A9) 


A.1.3  Tail  Width 

For  the  Fermi  function,  the  diffusivity  b  measures  how  fast  the  pd  curve  decays  to  zero  at 
ranges  beyond  R0.  Another  measure  of  this  that  is  generalisable  to  other  functions  is  the 
'width  of  the  tail  region'  Wt.  Many  definitions  are  possible,  but  for  present  purposes  it  is 
convenient  to  define  Wt  as  the  distance  from  R()  —  the  range  at  which  pd  =  0.5  — to  the  range 
at  which  pd  =  0.1,  as  illustrated  in  Figure  Al. 

For  the  Fermi  function,  substitution  in  Equation  (A8)  gives 


Figure  Al:  Definition  of  the  width  Wf  of  the  tail  region  for  a  pd  curve 
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log 


Wt  =Rnb- 


10po-l 

2p0-l 


0  l  +  blog(2p0-l)' 
For  exponentials,  substitution  in  Equation  (A4)  gives 


Wt  =  R( 


1-log2 

l°g(2Po)  ’ 


(A10) 


(All) 


Setting  Po  =  1  simplifies  both  these  equations  very  considerably.  In  the  case  of  the  Fermi 
function,  one  obtains 

Wt  =  R0&log9  =  0.9542R0&.  (A12) 

That  is,  Wt  is  proportional  to  b.  For  exponentials,  a  single  value  of  course  results: 

Wt  =  R0^^  =  2.322  R0.  (A13) 

log  2 


A.2.  Some  Results  from  Geometry 

A.2.1  Cassianian  Ovals 

A  'Cassinian  oval'  is  the  locus  of  a  point  moving  such  that  the  product  of  its  distances  to 
two  fixed  points  remains  constant.  That  is,  in  terms  of  the  quantities  shown  in  Figure 
A2(a),  the  product  Rst^TR  is  constant.  In  labelling  the  examples  in  Figure  A2,  we  set 

%rRTR=R02'  (A14) 

corresponding  to  the  50%  probability  contour  in  a  coverage  plot. 

Compilations  of  the  properties  of  Cassinian  ovals  are  available  (e.g.  [33(§5.16)]).  Their  be¬ 
haviour  divides  into  the  three  classes  shown  in  Figure  A2.  For  all  three,  the  maximum 


Figure  A2:  Examples  of  the  three  classes  of  Cassinian  oval:  (a)  RSK  <  R0V2,  (b)  R0V2  <  Rs^  <  2R0, 
(c)  Rsr  >  2R0.  The  'centre'  of  the  oval  is  the  midpoint  of  the  line  joining  the  foci,  even  if  this  does 
not  lie  within  a  loop,  as  in  panel  (c). 
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extent  xmax  of  the  oval  along  the  line  through  the  foci,  measured  from  the  centre,  is 


•^max 


(A15) 


The  classification  of  an  oval  depends  on  the  value  of  RgR  compared  with  R0 : 

•  Rsr  <  R0a/2  (Fig.  A2a):  The  oval  is  similar  in  appearance  to  an  ellipse.  The  maximum 
extent  ymax  from  the  line  joining  the  foci  occurs  at  the  perpendicular  through  the 
centre.  Its  value  is 

ymax=VRo2-RSR2/4-  (A16) 

•  RqV2  <  Rsr  <  2R0,  (Fig.  A2b):  The  oval  remains  in  one  loop,  but  the  maximum  ex¬ 
tent  ymax  now  occurs  away  from  the  perpendicular  through  the  centre.  Its  value  is 

y max  =  R02  / RSR  •  (  A1 7) 


•  Rg^  >  2  R(),  (Fig.  A2c):  The  oval  has  two  loops,  but  the  value  of  ymax  is  still  given  by 
Equation  (A17).  The  loops  intersect  the  line  joining  the  foci  at  ±  i'mm,  which  is 


•^min 


(A18) 


A.2.2  Ellipses 

The  ellipse  is  a  very  well  studied  curve,  but  its  properties  are  usually  given  in  terms  of  the 
maximum  extents,  commonly  called  the  major  and  minor  axes  — the  quantities  ,rmax  and 
3/max  Figure  A3  (e.g.  [33(§3.4)]).  For  present  purposes,  we  need  quantities  expressed  in 
terms  of  RgR  and  the  monostatic  blind-zone  radius  Rg,  so  as  to  understand  the  behaviour 
of  the  ellipse  as  these  two  parameters  change. 

As  is  well  known,  an  ellipse  is  the  locus  of  a  point  moving  such  that  the  sum  of  its 
distances  to  two  fixed  points  remains  constant.  To  represent  the  blind  zone,  we  write  (cf. 
Eqn  8) 

Rst  +  Rtr  =  Rsr  +  2Rb.  (A19) 

Standard  expressions  for  the  properties  of  the  ellipse  can  be  rearranged  in  favour  of  Rsr 
and  Rb  to  give 

-Lnax  -  +  ^Sr/^'  */max=  +  ^5R/^b  ■  (A20) 
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The  interesting  aspects  for  understanding  the  behaviour  of  the  blind  zone  are  the  beha¬ 
viour  as 

•  Rsr  increases,  for  constant  and 

•  decreases,  for  constant  Rsr; 

for  the  first  describes  the  growth  of  the  blind  zone  as  source  and  receiver  are  separated 
and  the  second  describes  the  effects  of  increasing  pulse  compression  at  fixed  source- 
receiver  separation.  These  are  depicted  in  Figure  A4. 


For  reference,  the  eccentricity  e  of  the  ellipse  is 

e-  RSR 
rSR  +  2Rb 


(A21) 


An  ellipse  becomes  more  circular  as  e  — >  0.  This  corresponds  to  Rsr  — >  0,  i.e.  monostatic 
geometry.  As  R^  — s >  0,  e  — ►  1,  which  means  that  the  ellipse  becomes  more  elongated;  it  col¬ 
lapses  toward  the  line  segment  joining  source  and  receiver,  as  illustrated  in  Figure  A4(b). 


A.3.  Field  Track-Initiation  Probability 

A.3.1  Listing  of  Formulae 

The  full  derivation  of  formulae  for  the  field  track-probability  pti  f  for  the  various  network 
architectures  with  the  3-in-5  track-initiation  rule  is  detailed  in  a  companion  report  [34]. 
Here  the  results  are  listed  for  convenience. 

Monostatic,  Distributed  Tracking 

Pttf  =i-n(i-Pty)/  (A22) 

7=1 

where  pti  ,•  is  the  track-initiation  probability  of  sensor  j,  given  by  Equation  (9),  and  J  is  the 
total  number  of  sensors  in  the  field. 


Figure  A4:  Blind-zone  behaviour  for  (a)  constant  =  0.2  and  source-receiver  separations  RSK 
shown  and  (b)  fixed  source-receiver  separation  Rs^  =  2.0  and  Rj,  values  shown.  Source  location  is 
marked  with  a  star  and  receiver  location  with  a  circle. 
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Monostatic,  Centralised  Tracking 


The  field  track-initiation  probability  is  written  as: 

pti,{  =  l-p(0)-p(l)-Pl(2)-p2(2),  (A23) 

where,  after  five  pings  of  each  source  in  the  field,  p(0)  is  the  probability  that  there  have 
been  no  detections  anywhere,  p(  1)  is  the  probability  of  exactly  one  detection  somewhere  in 
the  field,  p1(2)  is  the  probability  of  exactly  two  detections  by  the  same  receiver  and  p2( 2)  is 
the  probability  of  exactly  two  detections,  but  by  different  receivers.  The  expressions  are: 

P(0)  =  n(l-Pd,,)5,  (A24) 


where  pd  j  is  the  detection  probability  for  sensor  ;  and  J  is  the  total  number  of  sensors  in 
the  field, 

p(l)  =  5p(0)f-^|-,  (A25) 


and 


J  (  ^ 

n( 2)=ioP(0)X 

jM-Pd.j) 


P2<?)  =  2 5P(0)X  t  7 - - wr - r 

j=l  M+t  l1  Vd Pd,j') 


Pd,jPd,j> 


(A26) 


(A27) 


Multistatic,  Distributed  Tracking 


As  with  the  monostatic  distributed  case,  each  receiver  is  first  treated  separately.  The  track- 
initiation  probability  for  receiver  k  is 

Pti,k  =  1  -  Pk( 0)  -  PicO)  -  Pt,k(2)  -  P2,k(2)-  (A28) 

The  probabilities  of  zero  and  exactly  one  detection  after  all  J  sources  in  the  field  have 
pinged  five  times  each  are 

w<o)=n(i-M5.  <A29> 

(A30) 

j=-\  1  ~  Pd,jk 


where;  enumerates  sources  and  p ^  is  the  detection  probability  for  source;  and  receiver  k. 
The  quantity  /q  ^(2),  for  which  both  detections  are  made  on  pings  from  the  same  source,  is: 


Pl,k(2)  =  10pfc(0) 


T  (  \ 

f  Pd,jk 
j=l[1~Pd,jky 


(A31) 
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and  the  expression  for  p2  /f(2),  for  which  the  two  detections  involve  different  sources,  is: 


P2,k(2)  =  25Pl(°)T.  L  77- 

HM+n 


Pd,jk  Pd,j'k 
Pd,jk)(^-Pd,j'k) 


(A32) 


Finally,  the  field  track-initiation  probability  ptif  is  obtained  from  a  combination  of  the 
track-initiation  probabilities  for  each  receiver: 

pti,f=i-n(i-ptu)-  <A33) 

k=l 


Multistatic,  Centralised  Tracking 


In  this  case,  pti  f  is  written  as 

Pti,f  =  1  -  P( 0)  -  P(l)  -  Pi(2)  -p2(2)  -  p3(2)  -  p4(2),  (A34) 

where  p(0)  and  p{  1)  are  the  probabilities  of  exactly  zero  and  one  detection  respectively  in 
five  pings  of  all  sources,  as  before,  p4(2)  is  the  probability  that  a  given  source-receiver  pair 
makes  exactly  two  detections,  p2(2)  is  the  probability  that  the  two  detections  involve  the 
same  source  but  different  receivers,  p3(2)  the  same  receiver  but  different  sources  and,  for 
p4(2),  both  sources  and  receivers  are  different  for  the  two  detections.  The  expressions  for 
the  probabilities  are: 

p<°)=nrD-M5,  <a35> 

j=lk=l 


where  is  the  detection  probability  for  the  pair  of  source  j  and  receiver  k,  which  may  or 
may  not  be  collocated, 

Ki)=5P(o)ii 


p,(2)  =  10p(0)f  f 

j= 1  k= 1 


and 


J  K-l  K 

p2(2)  =  25p(0)£  I  I  - 
j= 1  fc=lfc'=fc+l( 

7-1  /  fc 

p3(2)  =  25p(0)X  X  Xi 
7=1  j'=j+lk=l  { 

7  7  K-l  K 

p4(2)  =  25p(0)£  I  I  I 
7=1  /=i,  fc=i  L=fc+ 
(M) 


Pd,jk 

(A36) 

-i —  Pdq'fc 

/  X 

Pd,jk 

2 

(A3  7) 

y1~Pd,jk^ 

/ 

Pd,jk  Pd,jk' 

(A38) 

-i  —  Pd,jk)(^  ~  Pd,jk'^) 

Pd,jk  Pd,fk 

(A39) 

1  Pd,jk){l 

Pd,j'k) 

Pd,jk  Pd,j'k' 

(A40) 

1  (l  ~  Pd,jk)(^  "  Pd,j'k')  ' 

In  Equations  (A36)-(A40),  p(0)  refers  to  Equation  (A35),  not  Equation  (A24).  For  the  field 
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layout  used  in  this  report,  p2(2)  =  ^3(2)  because  sources  and  receivers  are  collocated,  but 
this  would  not  otherwise  generally  be  the  case. 

A.3.2  Restriction  to  Five  Pings 

The  track-initiation  probability  given  by  Equations  (A22),  (A23),  (A28)  or  (A34),  as  the  case 
may  be,  is  the  probability  of  starting  a  track  after  exactly  5  pings.  Usually  many  more 
pings  than  this  would  be  emitted  during  a  tactical  scenario.  This  can  be  handled  analy¬ 
tically  by  a  sliding-window  approach,  such  as  that  discussed  in  §4.3  of  reference  12  and 
used  in  the  simulation  described  in  reference  26.  However,  sliding-window  methods  are 
not  appropriate  for  the  present  study  because  there  is  no  natural  upper  limit  to  the  num¬ 
ber  of  pings  that  should  be  considered,  and  all  probabilities  approach  1.0  as  the  number  of 
pings  increases  [12],  The  impact  of  this  limitation  on  the  number  of  pings  considered  when 
computing  the  measure  of  performance  is  mitigated  by  the  fact  that  the  limitation  applies 
equally  to  all  cases. 
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Appendix  B: 

Elaborations  of  the  Method 

B.l.  Behaviour  at  Large  Source-Receiver  Separation 

Section  2.2.1  illustrates  a  correction  applied  to  the  definition  of  Rec|uiv  to  overcome  a 
problem  with  the  behaviour  of  Equation  (5)  when  the  target  is  close  to  a  receiver  that  is 
itself  distant  from  the  source.  As  discussed  in  §1.2.3,  there  are  three  candidates  for  the  cor¬ 
rection:  a  variation  in  the  spreading  law,  near-field  effects  and  absorption.  The  first  two  do 
not  work,  for  the  reasons  given  in  §§B.1.3  and  B.1.4  respectively.  Subsections  B.1.1  and 
B.1.2  detail  the  one  that  does  work:  absorption. 

B.1.1  Correcting  for  Absorption 

With  the  inclusion  of  absorption.  Equations  (3)  and  (4)  for  transmission  loss  become 

TLmono  =  201°SK2  +  2aK 

TLmulti  =  20  log  (RST  RTr)  +  a(RST  +  RTr)'  (B1) 

where  a  is  the  attenuation  coefficient.  It  is  clearly  not  possible  to  formulate  a  single  defi¬ 
nition  of  Requiv  that  covers  both  the  spreading  and  absorption  terms  simultaneously.  We 
expect  that,  at  acoustic  frequencies  typical  of  ASW  sonar,  spreading  will  dominate  except 
at  large  distance,  so  we  retain  Equation  (5)  as  the  definition  of  Ret]uiv  and  apply  a 
correction  for  absorption,  using  the  so-called  'transition  curve'  [29(Fig.l2.10)]  to  relate 
transmission  loss  and  detection  probability.  The  calculation  involves  the  following  steps: 

•  The  uncorrected  multistatic  signal  excess  SEuncorr  is  computed  from  the  uncor¬ 
rected  multistatic  detection  probability  Rd(^equiv)  hy 

SEuncorr  =  unV2  inverf  (2pd  ( 

^equiv )  ^  ),  (B2) 

where  inverf(-)  is  the  inverse  of  the  error  function  (eg.  [45(ch.7)])  and  an  is  the  stan¬ 
dard  deviation  of  the  intensity  of  oceanic  ambient  noise.  Following  Urick  [29],  we 
take  an  =  8  dB. 

•  Since  transmission  loss  contributes  linearly  (in  dB  units)  to  signal  excess,  allowance 
for  absorption  amounts  to  a  linear  correction  involving  the  difference  in  path 
lengths  in  the  two  lines  of  Equation  (Bl): 

SEcorr  —  SEuncorr  —  a ( Rgq  +  %R  —  2^equiv)  •  (h3) 

The  last  term  recognises  that  the  monostatic  pd  value  includes  a  contribution  from 
absorption,  which  means  that  the  condition  for  the  correction  to  be  zero  is  stronger 
than  Rsr  =  0,  it  encompasses  all  situations  where  the  target  is  equidistant  from  the 
source  and  receiver  (i.e.  RST  =  Rtr).  The  range  difference  —  the  term  in  brackets  in 
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Equation  (B3)  — can  never  be  negative  (Appendix  B.1.2);  that  is,  SEcorr  <  SEuncorr 
always. 

•  Finally,  the  corrected  signal  excess  is  converted  back  to  a  detection  probability: 


Pd,corr 


1 

f  SE  Y 

1  +  erf 

^^corr 

2 

v 

(B4) 


B.1.2  Notes  on  the  Absorption  Correction 
Conversion  to  and  from  Signal  Excess 

The  relationship  between  detection  probability  and  signal  excess  (Eqns  B2,  B4)  adopts  the 
standard  interpretation  that  SE  =  0  is  equivalent  to  pL\  =  0.5,  derived  from  the  definition  of 
signal  excess  as  the  amount  by  which  the  received  signal  level  exceeds  an  operator-speci¬ 
fied  detection  threshold.  This  is  in  accordance  with  the  standard  conception  of  a  received 
signal  as  composed  of  a  deterministic  part  arising  from  acoustic  scattering  from  external 
objects  and  a  randomly  fluctuating  component  arising  from  ambient  oceanic  noise  [35]. 
Then,  given  the  reasonable  assumption  that  the  random  fluctuations  are  symmetric  about 
their  mean,  one  expects  that,  when  the  average  signal  level  exactly  equals  the  detection 
threshold,  the  total  signal  will  exceed  the  threshold  for  half  the  time  and  fall  below  it  for 
the  other  half.  Hence  pd  =  0.5. 

Obtaining  SE  values  at  pd  values  away  from  0.5  requires  further  assumptions  on  the  nature 
of  the  oceanic  noise  fluctuations.  Equations  (B2)  and  (B4)  assume  Gaussian  noise  with 
standard  deviation  of  an.  There  is  less  empirical  support  for  this  assumption  than  for 
symmetry  about  the  mean.  However,  as  elsewhere,  we  rely  on  the  relative  nature  of  the 
analysis:  the  assumptions  impact  the  results  only  to  the  extent  that  they  affect  the  mono¬ 
static  and  multistatic  geometries  differently. 


Correction  to  Signal  Excess  is  never  Positive 

We  seek  to  demonstrate  that,  in  Equation  (B3),  SEcorr  is  not  greater  than  SEuncorr,  irrespec¬ 
tive  of  the  geometry.  It  amounts  to  showing  that 

D  =  Rst  +  Rtr  -  2Requiv  >  0  (5) 

always.  Define  the  following  quantities: 

R  =  2(kst  +  rtr),  a  =  2Irst~rtrI.  (6) 

Then,  because  of  its  symmetry,  the  expression  for  D  is 

D  =  2R-2Ryl(l-A)(l  +  A)/  (7) 

regardless  of  which  of  RST  or  Rtr  is  the  larger.  Hence 

D  =  2R(1-Vi-A2)/  (8) 
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which  is  not  less  than  zero  for  any  value  of  A.  This  is  hardly  surprising,  for  the  analysis 
above  merely  rehearses  the  proof  that  the  geometric  mean  of  two  real  numbers  is  never 
larger  than  their  arithmetic  mean— a  very  well  known  result  (e.g.  [46(p.504)]). 

Wlien  pd  Equals  1.0 

Equation  (B2)  gives  SEuncorr  =  co  when  pd  =  1.0.  It  means  that  the  correction  procedure 
does  not  work  with  an  ideal  definite-range  pd  curve.  Here  is  another  reason  for  going 
beyond  the  cookie  cutter  for  modelling  multistatics. 


B.1.3  Variation  in  the  Spreading  Law 


Equation  (5)  is  based  on  spherical  spreading,  which  assumes  a  uniform  ocean  in  all  direc¬ 
tions.  All  environments,  however,  are  limited  in  vertical  extent,  by  the  surface  and  by 
either  the  bottom  or  the  maximum  depth  of  acoustic  rays  in  upwardly  refracting  condi¬ 
tions.  The  effect  of  this  on  spreading  loss  is  usually  modelled  as  spherical  spreading  out  to 
a  'transition  range'  Rc,  beyond  which  the  spreading  is  cylindrical  [29(p.l52),47],  That  is. 
Equation  (3)  becomes 


|401ogR  if  R<RC 

}201og(RcR)  if  R  >  Rc 


(B9) 


and  its  multistatic  equivalent.  Equation  (4),  becomes 


TL 


multi 


20  log  (RstRtr) 

if  R^t  —  Rq  and  R'yr  —  Rq 

10  log  (RsTRtr2Rc) 

if  Rc^  >  Rq  and  R~[j^  —  Rq 

lOlog  (rst2RtrRc) 

if  Rst  —  Rq  ^^id  R'y r  >  Rq 

lOlog  ^RSTRTRRc2j 

if  Rc,t  ^  Rq  and  R^ r  >  Rq . 

(BIO) 


It  is  immediately  clear  that  cylindrical  spreading  cannot  provide  a  solution  to  the  dif¬ 
ficulties  with  Equation  (5);  for  all  the  lines  in  Equation  (BIO)  involve  the  product  of  RST 
and  Rtr,  so  all  exhibit  the  same  undesirable  behaviour  as  Equation  (5).  Nevertheless,  we 
added  Equation  (BIO)  to  the  MATLAB  code  to  check  its  effect,  by  adapting  the  procedure  of 
Appendix  B.  1.1.  In  summary: 

•  Compute  multistatic  pd  using  Equation  (5)  and  the  monostatic  pd  curve  of  choice. 

•  Convert  to  uncorrected  signal  excess  using  Equation  (B2). 

•  Apply  the  correction  as  the  difference  between  Equations  (B9)  and  (BIO). 

•  Convert  back  to  pd  using  Equation  (B4). 

Results  are  shown  in  Figure  B1  for  the  somewhat  extreme  case  of  a  large  source-receiver 
distance  — 9R0  — and  a  small  transition  range:  just  R0.  Despite  this  effort  to  magnify  the 
effect,  it  remains  modest,  as  indicated  by  dotted  line. 
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Figure  Bl:  As  in  Figure  6,  but  with  a  dotted  blue  line  showing  the  effect  of  cylindrical  spreading 
assuming  Rc  =  R0  (and  no  absorption  correction) 

Note  also  that  this  correction  goes  in  the  opposite  direction  to  the  correction  for  absorp¬ 
tion;  that  is,  it  increases  coverage  area.  This  is  to  be  expected;  for,  with  cylindrical  spread¬ 
ing,  signal  level  decreases  more  slowly  as  range  increases  than  with  spherical  spreading. 

B.1.4  Near-Field  Effects 

Expressions  such  as  Equations  (3),  (4),  (B9)  and  (BIO)  are  accurate  only  beyond  a  certain 
distance  from  a  scatterer,  in  what  is  known  as  the  Tar  field'.  Far-field  conditions  involve: 

♦  the  range, 

♦  the  wavelength  X  of  the  acoustic  wave, 

♦  the  shape  of  the  scatterer  (i.e.  sphere,  cylinder,  etc.),  and 

♦  the  dimensions  and  orientation  of  the  scatterer. 

The  condition  on  wavelength  is  X  <  2 na,  where  a  is  a  'typical'  dimension  of  the  scatterer 
[29(p.303)].  At  a  frequency  of  2  kHz,  the  acoustic  wavelength  in  water  is  about  0.8  m,  and 
smaller  at  higher  frequencies,  hence  submarine-sized  objects  fulfil  this  condition. 

The  condition  involving  range  depends  on  the  shape  and  orientation  of  the  scatterer.  For  a 
sphere,  one  requires  R  >  a,  where  a  is  the  radius  of  the  sphere  [29(p.303)j.  This  also  applies 
to  a  part  of  a  sphere.  So,  for  example,  the  bow  section  of  a  submarine  may  be  piecewise 
spherical  with  a  radius  of  several  metre.  Then  far-field  conditions  apply  that  same  number 
of  metres  from  the  bow.  In  other  words,  when  viewing  the  submarine  bow-on,  the  near¬ 
field  region  is  negligibly  small,  from  the  point  of  view  of  ASW  operations. 

On  the  other  hand,  the  far-field  condition  perpendicular  to  the  axis  of  a  cylinder  of  length 
L  and  radius  a  is  R  >  L2/X  (and  also  X  <  2 na,  as  above).  For  a  length  of  80  m  and  X  =  0.8  m 
(i.e.  an  acoustic  frequency  of  2  kHz),  this  evaluates  to  R  >  8  km,  indicating  that  near-field 
effects  may  be  present  at  tactically  interesting  ranges  when  beam  on  to  a  submarine. 
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The  near  field  can  be  thought  of  as  a  region  of  overlapping  wavelets  that  spread  spheri¬ 
cally  from  every  point  on  the  surface  of  the  scatterer.  It  means  that  the  acoustic  field  may 
display  interference  minima  and  maxima,  the  amplitude  and  locations  of  which  cannot  be 
determined  simply.  The  most  that  simple  arguments  can  give  are  the  expectations  that  the 
amplitude  should  diminish  with  increasing  distance  from  the  scatterer,  to  be  negligible  at 
the  range  given  by  the  far-field  condition,  and  that  interference  features  should  be  spaced 
about  X  apart,  which  is  of  order  meters.  This  spacing  makes  it  difficult  to  fall  back  on  our 
usual  argument  that,  because  the  effect  applies  equally  to  monostatics  as  multistatics,  it 
ought  to  largely  cancel;  for  moving  a  receiver  by  a  distance  only  of  order  X  may  change  the 
detection  probability,  perhaps  markedly,  depending  on  the  interference  amplitude  in  the 
region  in  question.  However,  this  picture  also  suggests  that  near-field  effects  are  (approx¬ 
imately)  periodic,  which  means  that  they  (approximately)  average  out  over  distances  of 
order  many  X,  which  in  turn  would  suggest  that  they  ought  not  to  have  a  large  impact  on  a 
total  coverage  area.  Real  contours  in  plots  like  Figure  11  may  be  less  regular  than  shown, 
leading  to  small  (dimensions  less  than  or  of  order  A.)  undetected  coverage  gaps,  but  the 
effects  probably  wash  out  when  averaged  over  regions  of  dimension  much  larger  than  X, 
and  would  in  any  case  be  broadly  similar  for  monostatics  and  multistatics.  Hence  we 
expect  any  effects  on  the  final  comparative  results  to  be  small. 

However  much  the  above  argument  may  seem  unsatisfactory  (because  unquantitative),  it 
is  the  best  currently  available,  to  our  knowledge,  short  of  a  full  coherent  wavelet  summa¬ 
tion,  and  that  would  take  the  calculations  to  a  level  of  complexity  scarcely  compatible  with 
the  aspirations  of  operations  analysis. 


B.2.  Bearing-Dependent  and  Realistic  Detection  Probabilities 

The  monostatic  detection-probability  models  used  in  this  report— Equations  (1),  (6)  and  (7) 
—  are  schematic.  It  is  clearly  straightforward  to  extend  the  method  to  realistic  pL\  curves, 
such  as  those  from  acoustic-propagation  modelling  (e.g.  [12(§4.6)]).  This  may  provide  a 
method  of  approximating  a  full  multistatic  acoustic-propagation  calculation,  which  can  be 
very  time-consuming.  However,  the  degree  of  accuracy  achievable  is  yet  to  be  explored 
and  is  probably  environment-dependent. 

Equations  (1),  (6)  and  (7)  are  also  independent  of  the  bearing  of  the  target  from  the  sensor. 
This  is  appropriate  for  sonobuoys  and  dipping  sonars,  but  towed  line  arrays,  for  example, 
show  lower  pd  at  end-fire  than  broadside.  It  is  clearly  not  difficult  to  extend  the  analysis  to 
bearing-dependent  sensors.  The  coverage-area  results  would  then  depend  on  the  relative 
orientation  of  neighbouring  sensors  as  well  as  the  spacing  between  them.  An  example  for  a 
small  field  of  monostatic  sensors  is  shown  in  Figure  33  of  reference  12. 

B.3.  Aspect-Dependent  Target  Strength 

The  possibility  of  'glint'  (i.e.  specular)  reflections  from  a  target  is  one  of  the  attractions  of 
multistatics,  and  there  are  analyses  of  how  to  lay  out  a  field  to  maximise  it  (e.g.  [22]).  Glint 
reflections  can  occur  with  monostatic  sensors,  of  course,  but  submarine  COs  are  careful  to 
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avoid  them,  which  they  can  do  because  they  know  the  bearing  of  each  active  sensor  within 
range.  When  facing  a  multistatic  field,  on  the  other  hand,  the  submariner  knows  only  the 
locations  of  the  sources. 

The  method  of  correcting  for  absorption  provides  a  template  for  the  inclusion  of  aspect- 
dependent  target  strength:  having  converted  pL\  to  signal  excess,  a  correction  for  bistatic 
target  strength  can  be  applied.  In  general,  bistatic  target  strength  is  a  function  of  two 
angles,  the  bistatic  angle  cp  and  the  bistatic  aspect  angle  a,  as  defined  in  Figure  2  (p.  4).  The 
so-called  'bistatic  theorem'  (e.g.  [48])  states  that  target  strength  is  a  function  of  a  only,  but 
this  is  an  approximation,  as  Cox  explains  [15]. 

Inclusion  of  aspect-dependent  target  strength  means  that,  on  a  coverage-area  plot,  the 
probability  becomes  a  function  of  a.  Hence,  some  assumption  about  submarine  location 
and  heading  is  required.  The  results  do  not  explicitly  depend  on  cp,  regardless  of  whether 
the  bistatic  theorem  is  adopted,  because  cp  has  a  unique  value  at  each  point  on  a  coverage- 
area  plot.  That  is,  such  cp  dependence  as  may  exist  and  as  the  analyst  chooses  to  include  is 
automatically  accounted  for  in  the  coverage  calculations. 


36 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-TR-2562 


Appendix  C: 
Matlab  Code 


This  appendix  contains  listings  of  two  Matlab  codes  used  to  compute  the  multistatic 
coverage-area  diagrams  in  the  body  of  this  report.  The  first,  giving  for  a  bistatic  pair  in¬ 
cluding  both  the  absorption  and  cylindrical  spreading-loss  corrections,  produced  Figures 
3,  5-8  and  Bl.  The  second,  giving  field  track-initiation  probability  for  a  field  of  collocated 
source-receiver  pairs,  includes  only  the  absorption  correction.  It  produced  Figures  11-15 
and  the  results  in  Figures  16  and  17. 


C.l.  Bistatic  Pair 

%  Plotting  &  display  code  for  the  bistatic  case  when  transmission  loss  is  dealt  with  properly  for 
%  absorption  AND  spreading,  showing  contours  of  DETECTION  PROBABILITY  (not  pti). 

% 

clear  all 

%  - Data  entry - 

% 

%  Set  the  vector  of  source-receiver  separations;  the  code  generates  one  plot  for  each 
%  separation.  A  monostatic  calculation  results  when  Rsr  is  zero.  Rsr  is  a  fraction  of  rO. 

Rsr  =  [0  1  2  3  4  5  6  7]; 

%  set  the  type  of  the  equivalent  monostatic  pd  profile 

%ptype  =  'cookie'; 

%ptype  =  'exp'; 
ptype  =  'Fermi'; 


%  set  parameters  of  pd-profile  functions  (the  third  is  used  by  Fermi  only) 

p0=1.0;  %  peak  pd 

rO  =  10.0;  %  characteristic  range  in  km 

b  =  0.5;  %  diffusivity  (fraction  of  rO)  for  Fermi  function 

%  Set  the  pulse  length  Rp  and  multistatic  signal-processing  pulse-compression  factor  pcf. 

%  The  target  is  in  the  blind  zone  if  Rst  +  Rrt  <  Rsr  +  Rp/pcf  (multistatics)  or  R  <  Rp/2  (monostatics). 
%  (Rrt  —  receiver-target  distance,  Rst  —  source-target  distance.) 

Rp  =  0;  %  in  km;  set  to  zero  to  switch  blast-zone  checking  off. 

pcf  =  1;  %  must  be  greater  than  or  equal  to  1  (=  1  for  no  compression) 

%  Set  attentuation  coefficient  to  be  used  in  the  absorption  correction  &  noise  standard  deviation 

attenuation  =  0.0;  %  dB/km;  set  to  zero  to  switch  correction  off 
noise_sd  =  8.0;  %  dB 

%  Set  whether  to  include  cylindrical  spreading  transmission  loss.  TL  =  'spreading'  means 
%  include.  Set  Rc,  the  range  at  which  to  switch  to  cylindrical  spreading. 
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%TL  =  'no  correction'; 

TL  =  'spreading'; 

Rc  =  10;  %  km 

%  Set  the  scale  for  the  contour  plots  (in  km).  (The  first  line  is  the  default;  switch  to  the 
%  second  line  for  explicit  control  of  the  scale.) 

plot_scale  =  max(rO*Rsr); 

%plot_scale  =  50; 

% - End  data  entry - 


%  Set  up  grid  of  target  locations.  Receiver  is  located  at  the  origin,  source  is  on  the  x  axis  at 
%  location  +Rsr(i).  However,  all  i  contour  plots  have  the  same  plot  limits,  determined  by  the 
%  largest  of  the  Rsr  values  (or  the  nominated  value). 

if  strcmp('exp',  ptype) 

range_of_the_day  =  rO  *  Iog10(2*p0); 
around  =  range_of_the_day  *  1.5  ; 
elseif  strcmp('Fermi',  ptype) 

range_of_the_day  =  r0*(1  +  (1/b)*log10(2*p0  - 1)); 
around  =  range_of_the_day  *  2.3; 
else 

range_of_the_day  =  rO; 
around  =  range_of_the_day  *  1 .05; 
end 

meshstep  =  (plot_scale  +  around*1 .7)/1000; 
targx  =  -around:meshstep:plot_scale  +  around*0.7; 
targy  =  -around:meshstep:around; 

[gridx,gridy]  =  meshgrid(targx,  targy); 

%  Calculate  receiver-target  distances 
Rrt  =  sqrt(gridx.A2  +  gridy.A2); 

%  Determine  number  of  source-receiver  distances  to  plot;  set  contour  levels;  set  up  the  line 
%  plots,  including  constructing  the  title  &  labelling  the  axes 

Rsrjen  =  size(rO*Rsr); 
nplot  =  Rsrjen(2); 
contourlevels  =  0.0:0.01:1.0; 
line  plot  =  figure;  orient  tall 
subplot(2,1,1);  hold  on 

plot  title  =  {'\fontsize{1  IjBistatic  {\itp}_d  (incl.  absorption  &  cylindrical  spreading)';  ... 

['Range  of  the  day  =  '  num2str(range_of_the_day,%1 1 .2g') '  km_ '];  ... 
['\fontsize{9}Pulse  length  for  blast  zone  =  '  num2str(Rp)  ... 

'  km,_  pulse  compression  factor  =  '  num2str(pcf)];  ... 

'4th  line';  '5th  line'; '  ';  '\fontsize{10}Transect  along  receiver-source  line'}; 
if  strcmp('spreading',  TL) 

plot_title{4}  =  ['attenuation  coeff  =  '  num2str(attenuation) '  dB/km,_  noise  std  dev  =  '  ... 
num2str(noise_sd) '  dB,  transition  radius  =  '  num2str(Rc) '  km']; 

else 

plot_title{4}  =  ['attenuation  coeff  =  '  num2str(attenuation) '  dB/km,_  noise  std  dev  =  '  ... 
num2str(noise_sd) '  dB']; 
end 
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if  strcmp('exp',  ptype) 

plot_title{5}  =  ['Equivalent  monostatic  {\itp}_d:  exponential,  {\itp}_0  =  '  num2str(p0)  ... 
{\itr}_0  =  '  num2str(r0)  'km']; 
elseif  strcmp('Fermi',  ptype) 

plot_title{5}  =  ['Equivalent  monostatic  {\itp}_d:  Fermi,  {\itp}_0  —  '  num2str(p0) ... 

',  {\itr}_0  =  '  num2str(r0) '  km,  diffusivity  =  '  num2str(b)]; 

else 

plot_title{5}  =  ['Equivalent  monostatic  {\itp}_d:  cookie-cutter,  {\itp}_0  =  '  ... 
num2str(p0) ',  {\itr}_0  =  '  num2str(r0) '  km  ]; 
end 

title(plottitle) 

axis  ([min(targx),  max(targx),  0.,  1.0]) 
xlabel('Distance  (km)');  ylabel('{\itp}_d') 
xjnd  =  find((targy  +  meshstep/2)  >=  0,  1); 
yjnd  =  find((targx  +  meshstep/2)  >=  0,  1); 
subplot(2,1 ,2);  hold  on 

title('Transect  perpendicular  to  receiver-source  line  at  location  of  receiver') 
axis  ([min(targy),  max(targy),  0,  1.0]) 
xlabel('Distance  (km)');  ylabel('{\itp}_d') 
lineC  =  line_colour; 

%  For  each  source-receiver  distance  ... 
for  i  =  1:nplot 

%  compute  the  Rst  values  and  equivalent  monostatic  distances 
Rst  =  sqrt((gridx  -  rO*Rsr(i)).A2  +  gridy.A2); 

Requiv  =  sqrt(Rst.*Rrt); 

%  compute  the  pd  values  (ignoring  the  blast  zone  to  begin  with) 
if  strcmp('exp',  ptype) 

pd  =  pO  *  10.A(-Requiv/r0); 
elseif  strcmp('Fermi',  ptype) 

pd  =  pO  *  (1  +  10.A((Requiv/r0  -  1)/b)).A(-1); 
else 

pd  =(Requiv  <=  rO); 
pd  =  pO  *  pd; 
end 

%  set  up  for  corrections  —  compute  the  signal  excess  from  pd 
se  =  (sqrt(2)*  noise_sd)  .*  erfinv(2.*pd  -  1); 

%  apply  the  absorption  correction 

se  =  se  -  attenuation  .*  (Rst  +  Rrt  -  2*Requiv); 

%  check  whether  to  apply  correction  for  cylindrical  spreading 
if  strcmp('spreading  ,TL) 

%  Case  1  —  Rrt  >  Rc,  Rst  <=  Rc;  compute  the  cylindrical  transmission  loss 
%  minus  the  spherical  transmission  loss 
spreading_1  =  -20.*log10(Rst)  -  10.*log10(Rrt.*Rc); 

%  Case  2  —  Rst  >  Rc,  Rrt  <=  Rc;  ditto 
spreading_2  =  -20.*log10(Rrt)  -  10.*log10(Rst.*Rc); 

%  Calculate  regions  for  which  we  assume  equivalent  monostatic  is  in 
%  cylindrical  and  spherical  phases  and  add  on  the  appropriate  monostatic 
%  transmission  loss, 
sphere  =  (Requiv  <=  Rc); 
cylinder  =  (Requiv  >  Rc); 

spreading_1  (sphere)  =  spreading_1  (sphere)  +  40.*log10(Requiv(sphere)); 
spreading_1  (cylinder)  =  spreading_1  (cylinder)  +  20.*log10(Requiv(cylinder)*Rc); 
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spreading_2(sphere)  =  spreading_2(sphere)  +  40.*log10(Requiv(sphere)); 
spreading_2(cylinder)  =  spreading_2(cylinder)  +  20.*log10(Requiv(cylinder)*Rc); 

%  Calculate  regions  for  which  spreading_1  and  spreading_2  do  not  occur. 

outer_1  =  (Rrt  <=  Rc  |  Rst  >  Rc); 

spreading_1(outer_1)  =  0; 

outer_2  =  (Rst  <=  Rc  |  Rrt  >  Rc); 

spreading_2(outer_2)  =  0; 

%  Apply  correction 

spreading  =  spreading_1  +  spreading_2; 
se  =  se  +  spreading; 
end 

%  Recalculate  pd  from  signal  excess 
pd  =  0.5*(1  +  erf(se  ./  (sqrt(2)  *  noise_sd))); 

%  Blind  zone  —  set  up  the  logical  mask  indicating  which  target  locations  are  inside  the  blind 
%  zone  &  use  it  to  zero  out  that  part  of  the  pd  array.  Note  that  the  pulse-compression  factor 
%  does  not  apply  for  the  monostatic  configuration.  (Actually,  set  pd  values  to  -0.1  inside  the 
%  blind  zone  so  that  blind-zone  delineation  works  properly  in  the  plots.) 


Rb  =  Rp; 
if  (Rsr(i)  >  0) 

Rb  =  Rp/pcf; 
end 

blast  =  (Rst  +  Rrt  <  rO*Rsr(i)  +  Rb); 
pd(blast)  =  -0.1; 

%  add  lines  to  the  2-d  plot 

figure(line _ plot) 

subplot(2, 1,1) 

%  plot  pd  along  (or  near  to)  the  x  axis 

plot(targx,  pd(x_ind,:),  'LineStyle', 'Color',  lineC(i, :)) 

subplot(2,1,2) 

%  plot  pd  along  (or  near  to)  the  y  axis,  i.e.  receiver  location 
plot(targy,  pd(:,y_ind),  'LineStyle', 'Color',  lineC(i, :)) 

%  Set  up  for  the  contour  plot, 
figure;  orient  landscape;  hold  on 

axis  equal;  axis  ([min(targx),  max(targx),  min(targy),  max(targy)]); 
h  =  gca;  set  (h,  'TickLength',  [0.05,  0.1]); 
xlabel('Distance  (km)');  ylabel('Distance  (km)') 
caxis([0.0  1.0]);  colormap(pdmapMPF);  colorbar('EastOutside') 

%  Adapt  the  title  and  apply  ... 
plot_title{6}  =  '\fontsize{6}  '; 

plot_title{7}  =  ['\fontsize{1 1}_  Source  --  receiver  distance  =  '  ... 

num2str(rO*Rsr(i)) '  km_ ']; 
title(plottitle) 

%  Produce  a  filled  contour  plot,  then  overlay  with  a  thin  contour  line  around  the  blind  zone. 
%  Plot  markers  at  the  locations  of  the  source  and  receiver. 

contourf(gridx,  gridy,  pd,  contourlevels,  'LineStyle',  'none') 
contour  (gridx,  gridy,  pd,  [-0.1  -0.1],  '-k',  'LineWidth',  0.2) 
plot(rO*Rsr(i),  0,  'kp',  'MarkerFaceColor',  'k',  'MarkerSize',  40) 
plot(0,  0,  'ko',  'MarkerFaceColor',  'k',  'MarkerSize',  27) 
hold  off 
end 
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C.2.  Field  of  Collocated  Sources  and  Receivers 

%  Plotting  &  display  code  for  the  case  of  many  collocated  source-receiver  pairs. 

%  Can  be  run  distributed  monostatic,  networked  monostatic,  distributed  multistatic  or  networked 
%  multistatic.  Metric  is  probability  of  track  initiation  using  3-in-5. 
clear  all 


%  - Data  entry - 

%  Set  parameters  of  pd-profile  functions.  (The  third  is  used  by  Fermi  only.) 

p0=1.0;  %  peak  pd 

rO  =  10.0;  %  characteristic  range 

b  =  0.5;  %  (r_0.1  -  r_0.5)  /  r0*log10((10*p0  -  1)  /  (2p0  -  1))  ...  fraction  of  range  of  day 

%  Set  a  vector  of  sonar  spacing  values.  (One  plot  is  produced  for  each  value  in  such 

%  a  way  that  all  plots  have  the  same  scale.)  Then  set  the  array  of  source-receiver 
%  locations  as  a  fraction  of  that  spacing  value. 

%  (One  row  —  x  val,  y  val  —  for  each  sonar  in  the  field.) 

%area  =  [60]; 

spacing  =  [8.2  8.3  8.4]; 
pattern  =  [0  0;  ... 


0.0  1. 

0; 

0.0  -1.0;... 

0.86603 

0.5; 

0.86603  -0.5; 

-0.86603 

0.5; 

-0.86603  -0.5; 

0.86603 

1.5; 

0.86603  -1.5; 

-0.86603 

1.5; 

-0.86603  -1.5; 

1.73205 

0.0; 

-1.73205  0.0; 

1.73205 

1.0; 

1.73205  -1.0; 

-1.73205 

1.0; 

-1.73205  -1.0; 

0.0  2.0; 

0.0  -2.0 

]; 


%  Choose  one  of  the  following  4  lines  to  set  the  mode  of  operation  of  the  sonar  field. 
%mode  =  'Monostatic  with  distributed  tracking'; 

%mode  =  'Monostatic  with  networked  tracking'; 

%mode  =  'Multistatic  with  distributed  tracking'; 
mode  =  'Multistatic  with  networked  tracking'; 

%  Choose  one  to  set  the  type  of  the  equivalent  monostatic  pd  profile. 


ptype  =  'exp'; 
%ptype  =  'Fermi'; 


%  Set  the  blast-zone  length  Rb  and  multistatic  pulse-compression  factor  pcf.  The  target  is  in 
%  the  blast  zone  if  R  <  Rb/2  for  monostatics  or  Rst  +  Rrt  <  Rsr  +  Rb/pcf  for  multistatics. 

%  (R  =  monostatic  range,  Rst  =  source-target  distance,  Rrt  =  receiver-target  distance, 

%  Rsr  =  source-receiver  distance.) 

Rb  =  1 .5;  %  Distance  units;  set  to  zero  to  switch  blast-zone  checking  off. 

pcf  =  1 ;  %  Must  be  >=  1 ;  set  to  1  to  turn  pulse  compression  off. 

%  Set  the  scale  for  the  plot  (in  distance  units).  (The  first  line  is  the  default; 

%  switch  to  the  second  line  for  explicit  control  of  the  scale.) 
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plot_scale  =  max(rO*spacing); 

%plot_scale  =  20; 

%absorption  =  'off'; 
absorption  =  'on'; 

%Set  attenuation  (dB/km)  to  be  used  in  absorption, 
attenuation  =  0.1; 

%Set  noise  standard  deviation  (dB)  to  be  used  in  absorption. 
noise_sd  =  8; 


% - End  of  data  entry - 

% 

% 

%  Set  up  grid  of  target  locations, 
meshstep  =  plot_scale/50; 
max_vals  =  max(pattern)  *  plot_scale; 
min_vals  =  min(pattern)  *  plot_scale; 
if  strcmp('exp',  ptype) 
range_of_the_day  =  rO; 
around  =  range_of_the_day  *  3; 
else 

ra ng e_of_th e_d ay  =  r0+(r0/b)*log10(2*p0  - 1); 
around  =  range_of_the_day  *  1.5; 
end 

targx  =  (min_vals(1)  -  around):meshstep:(max_vals(1)  +  around); 
targy  =  (min_vals(2)  -  around):meshstep:(max_vals(2)  +  around); 

[gridx,gridy]  =  meshgrid(targx,  targy); 

%  Determine  number  of  field  spacings  to  compute  &  number  of  sonars  in  the  field; 

%  set  the  levels  for  the  filled  contour  plots. 

spjen  =  size(rO*spacing); 

ncases  =  spjen(2); 

snrjen  =  size(pattern); 

nsnrs  =  snrjen(l); 

contourlevels  =  0.0:0.01:1.0; 

for  icase  =  1  meases 

%lnitialise  ratio  of  detection  store.  Each  cell(source,  receiver)  of  r_d_store  contains  the  r_d 
%  values  for  that  particular  source  and  receiver.  Within  each  cell  is  a  gridx,gridy  matrix. 
r_d_store  =  cell(nsnrs, nsnrs); 

%initialise  matrix  for  prob  of  getting  zero  detections  at  all  receivers  after  each  source  has 
%  pinged  5  times  —  multistatics 
prO  =  ones(size(gridx)); 

%initialise  matrix  for  getting  zero  detections  at  all  receivers  after  their  respective  sources 
%  have  pinged  5  times  —  monostatics 
pr0_mono  =  ones(size(gridx)); 

%initialise  pti  field 
pti_field  =  ones(size(gridx)); 

%The  above  must  reset  for  each  of  the  ncases 

%  Calculate  sonar  locations. 

sonars  =  pattern  *  rO  *  spacing(icase); 

for  irevr  =  1:nsnrs 

%  Compute  the  array  of  source-target  distances. 

Rrt  =  sqrt((gridx  -  sonars(ircvr,1)).A2  +  (gridy  -  sonars(ircvr,2)).A2); 


42 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-TR-2562 


%initialise  matrix  for  probability  of  0  detections  at  k  after  5  pings  from  each  source 
%Must  reset  for  each  sonar  as  receiver 
prO_k  =  ones(size(gridx)); 

if  strcmp('Monostatic  with  distributed  tracking',  mode)  ||  ... 

strcmp('Monostatic  with  networked  tracking',  mode) 

%calculate  equivalent  multistatic  range  (not  incl  TLabsorption  -  only  TLspreading) 
Requiv  =  sqrt((gridx  -  sonars(ircvr,1  )).A2  +  (gridy  -  sonars(ircvr,2)).A2); 

%  Compute  the  pd  values,  ignoring  the  blast  zone, 
if  strcmp('exp',  ptype) 

pd  =  pO  *  10  A(-log10(2)*Requiv/r0); 
else 

pd  =  pO  *  (1  +  10  A(((Requiv/rO)  -  1  )/b)).A(-1 ); 
end 

%  Compute  &  apply  the  mask  to  zero  out  probabilities  in  the  blast  zone, 
blast  =  (Requiv  <  Rb/2); 
pd(blast)  =  0.; 

if  strcmp('Monostatic  with  distributed  tracking',  mode) 

%Calculate  individual  pti  for  sonobuoy  then  increment  pti  field 
pti  =  pd  A3  .*  (10  -  15*pd  +  6*pd  A2); 
pti  field  =  pti  field  .*  (1  -  pti); 
end 

%increment  probability  of  no  detections  from  any  sonobuoy 
pr0_mono  =  pr0_mono.*(1-pd)  A5; 

%calculate  ratio  of  detection  to  no  detection  and  store  each  receivers. 
r_d  =  pd  ./  (1  -  pd); 
r_d_store(ircvr,  ircvr)  =  {r_d}; 

else 

for  isrce  =  1:nsnrs 

%  Compute  the  source-receiver,  receiver-target  and  equivalent  monostatic 
%  distances  (N.B.  the  monostatic  case  in  INCLUDED;  i.e.  a  sonar  can  act  both  as 
%  source  and  receiver.) 

Rsr  =  sqrt((sonars(isrce,1)  -  sonars(ircvr,1))A2  +  ... 

(sonars(isrce,2)  -  sonars(ircvr,2))A2); 

Rst  =  sqrt((gridx  -  sonars(isrce,  1)).A2  +  (gridy  -  sonars(isrce,2)).A2); 

Requiv  =  sqrt(Rst.*Rrt); 

%  Compute  the  pd  values,  ignoring  the  blast  zone, 
if  strcmp('exp',  ptype) 

pd  =  pO  *  1 0.A(-log1 0(2)*Requiv/r0); 
else 

pd  =  pO  *  (1  +  10.A(((Requiv/r0)  - 1  )/b)).A(-1 ); 
end 

if  strcmp('on',  absorption) 

%compute  the  signal  excess  for  pd 
se  =  sqrt(2)*noise_sd  .*  erfinv(2.*pd  -  1); 

%add  in  the  absorption  term 

se  =  se  -  attenuation. *(Rst  +  Rrt)  +  2*attenuation*Requiv; 

%recalculate  pd  from  signal  excess 
pd  =  0.5*(1  +  erf(se  ./  (sqrt(2)*noise_sd))); 
end 
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%  Compute  &  apply  mask  to  zero  out  probabilities  in  the  blast  zone, 
blast  =  (Rst  +  Rrt  <  Rsr  +  Rb/pcf); 
pd(blast)  =  0.; 

%increment  pr0_k  -  probability  of  no  detections  at  receiver  k  from  any  source  after 
%each  has  pinged  five  times. 
pr0_k  =  pr0_k.*(1  -  pd).A5; 

%ratio  of  detection  to  no  detection. 
r_d  =  pd  ./  (1  -  pd); 
r_d_store(ircvr,  isrce)  =  {r_d}; 
end 
end 

%increment  prO  -  the  probability  of  no  detections  at  any  receiver  after  each  source  has 
%  pinged  5  times. 
prO  =  pr0.*pr0_k; 

%This  is  where  the  multi  dist  stuff  belongs  -  calculate  the  pti  for  each  receiver: 
if  strcmp('Multistatic  with  distributed  tracking',  mode) 

%start  tracks  at 
r_d_1  =  zeros(size(gridx)); 
r_d_2_1  =  zeros(size(gridx)); 
r_d_2_2  =  zeros(size(gridx)); 

%rd  calculation 
for  u  =  1:nsnrs 

r_d_1  =  r__d_1  +  r_d_store{u,ircvr}(:,:); 
r_d_2_1  =  r_d_2_1  +  (r_d_store{u,ircvr}(:,:)).A2; 
end 

%Calculate  probability  of  one  detection 
prl  =  5*r_d_1  ,*prO_k; 

%calculation  for  probability  exactly  2  (1).  10  is  from  5C2. 
pr2_1  =  1 0*r_d_2_1  ,*prO_k; 

%calculation  for  probability  exactly  2  (2).  25  is  from  (5C1)A2. 
for  u  =  1:nsnrs-1 
for  w  =  u+1  :nsnrs 

r_d_2_2  =  r_d_2_2  +  r_d_store{u,ircvr}(:,:).*r_d_store{w,ircvr}(:,:); 
end 
end 

pr2_2  =  25*  r_d_2_2 .  *  p  r0_k ; 
pr2  =  pr2_1  +  pr2_2; 

pti  =  1  -  pr0_k  -  prl  -  pr2; 
ptifield  =  pti_field.*(1  -pti); 
end 
end 

if  strcmp('Multistatic  with  distributed  tracking',  mode)  ||  ... 

strcmp('Monostatic  with  distributed  tracking',  mode) 

pti_field  =  1  -pti  field; 
end 

%For  the  netowrked  multistatic  case 
if  strcmp('Multistatic  with  networked  tracking',  mode) 
r_d_1  =  zeros(size(gridx)); 
r_d_2_1  =  zeros(size(gridx)); 
r_d_2_2  =  zeros(size(gridx)); 
r_d_2_3  =  zeros(size(gridx)); 
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%rd  calculation  for  probability  of  detecting  exactly  one  and 
%exactly  two  from  one  pair 
for  u  =  1:nsnrs 
for  v  =  1:nsnrs 

r_d_1  =  r__d_1  +  r_d_store{u,v}(:,:); 
r_d_2_1  =  r_d_2_1  +  (r_d_store{u,v}(:,:)).A2; 
end 

for  v  =  1:nsnrs-1 
for  w  =  v+1:nsnrs 

r_d_2_2  =  r_d_2_2  +  r_d_store{u,v}(:,:).‘r_d_store{u,w}(:,:); 
for  x  =  1  :nsnrs 
if  x~=u 

r_d_2_3  =  r_d_2_3  +  r_d_store{u,v}(:,:).*r_d_store{x,w}(:,:); 
end 
end 
end 
end 
end 

%Calculate  probability  of  one  detection 
prl  =  5*r_d_1  ,*prO; 

%Calculate  probability  of  2  detections  from  same  source  and 
%receiver 

pr2_1  =  10*r_d_2_1; 

%calculation  for  probability  exactly  2  (2).  25  is  from  (5C1)A2. 
pr2_2  =  25*r_d_2_2; 

%calculation  for  probability  exactly  2  (3).  25  is  from  (5C1)A2. 
pr2_3  =  25*r_d_2_3; 

pr2  =  pr0.*(pr2_1  +  pr2_2  +  pr2_3); 

pti_field  =  1  -  prO  -  prl  -  pr2; 
end 

if  strcmp('Monostatic  with  networked  tracking',  mode) 
r_d_1  =  zeros(size(gridx)); 
r_d_2_1  =  zeros(size(gridx)); 
r_d_2_2  =  zeros(size(gridx)); 
for  u  =  1:nsnrs 

r_d_1  =  rdl  +  r_d_store{u,u}(:,:); 
r_d_2_1  =  r_d_2_1  +  (r_d_store{u,u}(:,:))  A2; 
end 

for  u  =  1:nsnrs-1 
for  v  =  u+1  :nsnrs 

r_d_2_2  =  r_d_2_2  +  r_d_store{u,u}(:,:).*r_d_store{v,v}(:,:); 
end 
end 

prl  =  5*r_d_1  .*prO_mono; 

pr2  =  (10*r_d_2_1  +  25*r_d_2_2).*pr0_mono; 

pti_field  =  1  -  prO_mono  -  prl  -  pr2; 
end 

%  Set  up  to  plot  the  field  pti. 

figure  ('PaperType',  'A3');  orient  tall;  hold  on 

axis  equal;  axis  ([min(targx),  max(targx),  min(targy),  max(targy)]) 

xlabel  ('Distance  (km)');  ylabel  ('Distance  (km)') 

caxis([0  1.0]);  colormap(pdmapMPF);  colorbar('EastOutside') 
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%  Construct  and  apply  the  appropriate  title  ... 

plot  title  =  {'\fontsize{12}  Collocated  sources  &  receivers  _ '2nd  line';  ... 

'  Contours  of  {\itp} _ tj  using  3-in-5  track  initiation';  ... 

['\fontsize{10}_  Pulse  length  for  blast  zone  =  '  num2str(Rb)  ... 

'  km,  pulse  compression  factor  =  '  num2str(pcf)];  '5th  line';  ... 

'\fontsize{3} ';  ['\fontsize{12}  Sonar  separation  =  ' ... 

num2str(r0  *  spacing(icase)) '  km,  Range  of  the  day  ({\itR}_0)=  '  ... 

num2str(range_of_the_day,  '%11.2g') '  km  '];  ... 

['Absorption  =  '  num2str(absorption)]}; 
plot_title{2}  =  num2str(mode); 

plot _title{5}  =  ['Equivalent  monostatic  {\itp}_d:'  num2str(ptype) ',  {\itp}' ... 

'_0  =  '  num2str(p0) ',  diffusivity  =  '  num2str(b)]; 

title(plottitle) 

%  Produce  a  filled  contour  plot,  then  overlay  with  a  thick  contour  line  at  the  0.5  level  and 

%  draw  markers  at  the  location  of  each  sonar 

contourf(gridx,  gridy,  pti  field,  contourlevels,  'LineStyle',  'none') 

contour  (gridx,  gridy,  pti_field,  [0.5  0.5],  '-k',  'LineWidth',  1.2) 

for  isnr  =  1:nsnrs 

plot(sonars(isnr,1),  sonars(isnr,2),  'ko',  'MarkerFaceColor',  'k',  'MarkerSize',  4) 
end 

hold  off 
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